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ABSIRACT 


The 9cattefcLng of phonons from a substittitioha^ impTukt^j" in 
cubic lattices is studied in detail with the help of Gieen function 
technique used by Lifshitz and Montroll end Potts. Exact solutions 
for sca'fetering amplitude eie found explicitly and the scattering 
amplitude is expanded in terms of partial waves characteilaed by the 
irreducibLe representations of the point group of the lattice. The 
scattering cross sections are found in the limit of long waves showing 
Rayleigh type of scattering in absence of resonances which are taken 
care of by the resonance denominators. 

Specific examples of the scattering of phonons from substitutional 
impurities in simple cubic, body centred cubic and face centred cubic 
lattices are worked out in detail. Curves are given showing the 
dependence of resonance frequencies on parameters characterizing 
changes in mass and force constant. The long wavelength behavior of 
the scattering process is discussed and the total cross sections are 
found using optical theorem in this limit. The dependence of the 
scattering cross sections on the direction of incidence is shown by 
noting the different partial wave contributions for different directions 
of incidence. The inverse of mean free path for phonons due to 
scattering by substitutional impurities is found in the long wavelength 
limit. It is shown to be proportional tofiO ^ with resonance 
possibilities and the interference terms between partial waves 
characterized by different irreducible representations are given 
explicitly. Phonon scattering from a substi'tlanal impurity in a 



diatomic eiaple cubic lattice is also discussed and cUrves are drawn 
showing the dependence of resonance frequencies upon the pateiaeters 
ch*arecterlzlng changes in mass and force constant. 

Ai the end, a model of body centred cubic lattice is considered 
in tdaich the polarization of phonons is taken into account in the 
description of the scattering process, which is shown to depend upon 
the direction of incidence as well as the state of" polariiatiori of 'the 
incident phonons. The dispersion law for the model does not go ovar 
to the Debye model in the limit of long waves and oonsecpaently the 
dependence of the number of partial waires contributing to the scattering 
cross sectim even in the limit of long waves depends upon the incident 
direction and polarization of incident phcsn^as. The result for the 
isotope defect case in the limit of long waves differs from that 
calculated by assuming Debye model by a numerical factor which depends 
ripon the direction of incidence, A comparlsoTj is made betvreen this 
and the so called scalar models, • h 
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CHAPTBR I 


introductio:j 

The problem of scattering of phonons bj a substitutional 
impuid. was first treated by I.M. Lifshitz^ ^ whose pioneering work^ 
on the problem of Dyn&mic Theory of Won^ideal crystal lattice is well 
known. It has subsequently been discussed by many workers'^^'^including 
the author • Most of the fund-amental concepts and methodological 

aspects of the problem are present in Lifshitz’s work. . It has been 
shown that the problem can be worked out exactly using Green function 
techniques. The perturbation due to impurity may be large but the range 
of this perturbation is assumed to be quite small within the crystal. The 
principal result is that the perturbed region surrounding the impurity 
appears as a separate system whose energy le-vels depend upon the energy 
spectrum of the host lattice. If the perturbation is sufficiently large, 
an inpurity level may lie outside the band, with a long lifetime since 
there are no states in the neighbourhood into which it can decay. On the 
other hand for moderate perturbations the impurity level will lie inside 
the band with a finite probability of decay. The latter level has 
®^f®cts on the scattering of phonons from the impurity'. 

As an exanple Llfahitz, considered the case of an isotope defect 
and the scattering of phonons from it. He obtained the expression for the 
hcuttering jampiitude by expanding the proper Green function for large 
distances. The behavior of this process in the limit of long waves are 
also discussed. Montroll arxi Potts ^ independently ' considered the 
problem in the same way as Lifshitz did. They also showed that the 
case of diatomic lattices can be handled by the formalism for mono-* 
atomic lattices by performing the so called transformatlon ^^i 
Id-fshliz noticed the fact that , by the introducticm of a substitutional 

« 111 « » 
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impurity the. translational symmetiy of the lattice is destroyed . and 
only the point group syBunetry remains, the idea of expanding the 
phonon scattering amplitude in terms of partial wares chaiacterised 
hy the irreducible repiesentations of the point group of the lattice 
was developed by others 

The rest of the work in this field has been mostly matters of 
detail. The literature is extensive bit the notations are uncorrelated, 
•sometimes to the extent of confusing the reader. Maradudin*^ has 
discussed the details and refinements of the problem on several 
occasions. Takenq^’®, Callaway®, Wagner^^,Kiumhansl^^^^’^ and Thoma and 
Ludwig have also considered the details of the problem with its 
implications on the evaluation of thermal conductivity due to scattering 
of phonons from substitutional impuritJLss . In many of these the idea 
of resonant scattering has been stressed. Krumhansl and Matthew^*^ 
considered the problem in a one dimensional model in a rather complete 
way, and have demonstrated interference between different partial 
waves in the expression for relaxation time. This has been elaborated 
for three dimensions in the present work. Elliot and Taylor^^ solved 
the problem by using double time Green function technique in an 
elegant way. Recently Klein given a completely gfeneral 

formulation of the problem within the framewoiic of the quantum 
mechanical formal theory of scattering. Several important aspects 
like thh optical theorem and orthonormality of the solutions have been 
discussed in this work, Wagner^^ has discussed the case of molecular 
impurities by introducing the molecular Green function. 

There has been considerable overlaping among the approaches 



and results of the vailous workers mentioned above. The overp.ll 
picture at present is that we have made substantial progress in 
understanding the details of the problem, but this by no means is 
complete and final. Experiments in this field has lagged behind theoiy 
and new experiments may reveal more details than has been known 
hitherto. The main source of experimental information on scattering 
of phonons from substitutional impurities at present is the low 
tenpereture thermal conductivity data. By analysing this. Walker and 
Pohl^^*^^ were able to infer the presence of resonance scattering 
leading to the deviation from Klemens' result^^. Baumann^^ and Klein^^ 
also arrived at similar conclusions. Very recently Huebener^^ has 
estimated the scattering cross section for phonons scattered from 
vacancies by measuring the change in thermoelectric power in Platinum. 
1 remark relevant to these experiments is that usually the impurities 
do not always go into substitutional sites. They sometimes go into 
interstitial positions so that the results have to be interpreted 
cautiously. 

There is a great deal of similarity among the problems 
■ of scattering of phonons, band electrons and spin waves from 
substitutional impurities. This has been explored by Callai^ay® in 
detail. Formalisms identical to that for the phonon scattering for 
the other two cases have been developed by various worke , 

In this dissertation scattering of phonons from substituti-. > 
onal impuidties in some three dimensional lattice models are conside- 
red in detail. The scattering amplitudes are expanded in terms of 
partial wayes characterised by different irreducible representations 
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of the point group of the lattice anfl the acattering oioss sections 
aio eiraluatefl in the liMt of long wares. The question of resonant 
scattering of phonons in the optical batjl of a dlatoinc simple cubic 
lattice is discussed. At the end a model of body centred cubic lattice 
is considered i„ which the polarisation of phonons is taken into 
account in the description of the scattering process. The question of 
ths dependence of the scattering process on the direction of incident 
Phonons is also erramlned. Finally the implications of a departure from 
Itebye model in -the limit of long wares is discussed. 



CHAPTER n 


GEKEML FORMJLATION 

2,1 Introduction . 

The theoiy of scattering of phonons from substitutional 

impurities started with efforts concentrated on isotopic impurities 

im-olving only a change' in mass. It was further believed that the first 

24 

Born aporoximation results were adequate to describe the process , 

but this assumption was soon aband.oned . The low tempera.tU3re 

thermal conductivity experiments by Klein^®, Walker and Pohl^^“^^ and 
25 

Bau:nfun" showed the inadequacy of first Bom approxiine.tion results 
and the necessity of including resonant scattering which can "be done 
by solving the problem exactly. They further showed that the rather ... 
simple case of isotope defect was also inadequate. One should consider 
the scattering of phonons from point defects characterized by changes 
in ma^ as well as force constant. Moreover one should look for exact 
solutions rather than perturbation approximations . We consider such 
an approach to the problem of scattering of phonons from a substitutional 
impurity in the present work. 

Some other investigators have considered this approach as 
already mentioned in the last, chapter. Ifeny other workers'^-^^^ have 
been intere-sted in the determination of only the localized modes due 
to a 'general substitutional impurity' by which we will mean, 
throughout this work, a substitutional impurity characterized by 
changes in mass as well as force constants. However this determination 
of localized modes can be achieved in a natural way from the solution 
of the general scattering problem. 
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' The simplest mcy’el for the scattering of phonons from e 
general substitutional impurity has been considered by Kiumhansl and 
Motthew^"^. Although they consider a one-dimensional problem, the 
qualitative conclusions d.e rived aie quite instructive and useful. 
Klein^'^ develops a t— matrix formalism which is applicable to three 
diiensional crystals. This has been extended further by his latest 
work^^ which contains many important theorems and an averaging 
procedure for a distribution of defects. It is a self consistent 
formalism applicable to many models and yiflds results that 
identical to those of the simple formalism described below ( when 
applied to the same models), except that it is perhaps easier to 
exbract numbers out of the latter than from the former formalism. 

Also the explicit detailed treatment of the resonance denominators 
in the latter formalism facilitates the investigation of resonances 
and localized modes. 

Before discussing the present formalism it is worthwhile 
to note that this scattering problem is analogous to the potential 
scattering in ordinary quantum mechanics except for two impoortent 
differences. It is a particular case of the more general problem of 
scattering of quasi-particles in ciystall±ne solids which has been 
studied by Callaway^. One of the differences is due to the fact that 
the dispersion law giving the wave number dependence of the single 
particle energy reveals the typical band structure involving both 
an upper and a lower bound to the energy. This permits the possibility 
of occuiswttceof bound states both above and below the optical bands 
and above the aconstic band, depending on the nature of the 
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scatterer. Moreover for the scattering process both the incident and 
the scattered particles have their adsiissible energies limited to the 
band, The wave number is also limited to lie within the first Brilloiin 
Zone# The other difference is related to the syrmetr:/ of the problem. 
The interaction due to the scatterer usually possesses the symmetry 
of the point group of the crystal rather than the spherical symmetry of 
the potentials studied in the theory of potential scattering. This 
feature necessitates a change in the method of partial wave analysis of 
the scattering amplitude which now has to be decomposed in terms of 
partial we've*' characterized by the irreducible representations of the 
point group of the lattice. 

The present formalism is perhaps mathematically simpler and 
less formal than those cited abrrve . The scattering amplitudes and 
their partial wave espensions are explicitly presented and tho optical 
theorem in the usual form is employed to demonstrate the non— central 
effects of a general perturbation. The solution is exact and. includes 
resonant scatterings. The long wavelength limit is found easily. The 
interference between terms corresponding to different partial waves 
can also be seen explicitly in the evaluation of relaxation times. 

^ 2»2 The Model 

The formalism developed here depends upon a number of 
assumptions. The first is the adiabatic approximation in which it is 
assumed that the electrons instantenously take up a configuration 
appropriate to that of the displaced nuclei and that the change of ■ 
energy of the electrons in the distorted lattice contributes to an 
effective inteiwuclear potential. This is quite plausible for the 
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vibrations in the aooustio moles that „e »ill deal vdth. The other 
important approximation is the harmonic approximation. Since the 
harmonic force constants assumed here may have a temperature dependence. 
»e are in effect assuming a pseudo-harmonic approximation which is 
adecfiate for low temperature properties. 

Here we ccnsider a monoatomic lattice and investigate oases 
where the question of polarization floes not enter into the picture. The 
30 called -scalar models- talong to this category. In addition there 
are models other than these scalar models which do not mix polarization 
such as the one discussed in the next chapter. The discussion of 
polarization is taken up in chapter y of the present work and a 
soussion of a diatomic latticB In the next chapter shows how the 
iBsults for the monoatomlo lattices can be easily carried over to the 
ease of diatomlo lattices by using the so called M* trensformation®. 

With slight modllioations, the results obtained here can -be carried 
over to the case of the scattering of spin waves ^ and the 
scattering of band electron ^ , 

e-Lectrons from substitutional impurities in 

the lattice . 

lattice under consideration has N atoms of mass M each 
which are ccuplsd to their nei^beurs threirgh harmonic force constants. 

The time Independent equations of motion can be written as, ' 

AH - (M.^)U = 0 ( 2.1 ) ^ 

Where A is the force constant matrix of dimensions N x N 
and 2 is the N-dlmensional vector whose elements are the one 
conponent displacements ( scalar field ) of atoms from their 
equilibrium sites lAth a time dependence of the fom exp(.^«jt ). 
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The matrix A can be shorn to be cyclic and hence its 
eigenvectors (k ) have elements of the form 

^0 (]£ »£ )= exp(lk«R ) (2.2 ) 


whore g denotes the lattice sites xrith leference to a given 
origin. The vector k ( with N different values) specifies the 
eigenvectors and the corresponding eigenvalues. Uq (k , R ) has the 
form of a running wave or a phonon propagating with a wave vector 
k * Substituting equation (2 '2) in .equation (2.1), one gets the 
eigenvalues of 2 which give the frequency squared as a function of 
k i,e the dispersion law which is written as u) (k ) and this has 
both an upper and a lower bound defining the band. The distinct 
values of k lie within the first Billlouin Zone. The eigenvectors 
S.O ( ) satisfy the orthogonality and closure properties in an 

obvious way 


and 


p, TIo*(k , S') Do (k , £ ) - pR , 2' 


23 S ) Do (k , H ) = ^k, k' 

R 


( 2.3 ) 


These properties of A and its eigenvectors are direct 
consequences of the cyclic boxondaiy conditions imposed on the 
lattice . 

. g.3.The perturbed lattice . 

Let a substitutional impurity be put into the lattice , This 
changes the mass of the atom at the substituted site to M + M 
and. the cbupling of this impurity atom to its nearest neighbours 
is given ly the changed harmonic force constant -j- C^Y 
/|\Ydenotes the change in the noimal value Y * Thus the impurity 
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affects the site it occupies plus its neighbours whose number is 
detennined by the range of interaction. JO.! these sites including 
that of the impurity atom are caUed the sites affected l^y the 
perturbation. In most of the cases these affected sites can be assumed 
to have the symmetiy of the point group of the lattice with the impurity 
si-oe ae the point left invariant under the opeiations of this point • 
group. Strictly speaking, the resulting symmetry group of the system 
will be the intersection of the point grrup of the lattice and the 
sjn^etxy group of the pertuxhation . What we are assuming here in 
Bct amounts to believing that the symmetry group nf the 

perturbation haa at laaat all the elamenta of the point group of 
the lattice. 

Then the displacements of all the atoms at the affected sites 
form a representation space of the point group of the lattice in Aich 
the reducible represertatlon of the Point group can le decomposed 
in terms of the irreducible representations of the point 

group as 

r "S (2.4) 

where are zero or integers. The above representation space 
will be called the -space. 

The equations of motion of the perturbed lattice can be 
written as 

42 -(Mu)2)D .^2 (2,5 ) 


*eie the K rt H perturbation matrte P has non- vanishing elenants 


only in ^ P.^aoe and can be partitioned into the form 
P.P O') 

°J ( 2-6 ) 
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wLaro the matrix p is of dimensions n x a and n is the number of 
iF-tticG sites affected by the perturbation^ i.e. the diuBnsionality of 
the p-^space. 

The solution of equation (2.5) outside the unperturbed band 

will give the discrete frequencies corresponding to bound states in 

potential scattering. They are the localised modes of the lattice 

with impurity. The scattering state solutions are obtained for value s‘ 

2 

of uj within the unperturbed band and can be written as satisfying 
an equation of the form , 

2 - Uq ^ (2.7 ) 

where the Green function matrix G is defined as , 

A ^ A P A —1 

G=-.(A-Ma)l) (2.8) 

where ? is the unit matrix. Uq is the unperturbed solution having 
the same eigenvalue as U and this is always possible in a large 
lattice^ because the frequencies form a continuum., 

The matrix elements of G can be written explicitly as^ 


0 (u? , E , ^)= Uo(k ) TJo(k,fi') 

" P / .1 \ 2 


60^(k)- 


(2.9 ) 


The summation over k is ever the first Billlouin Zone. 

The disadvantage of the use of such Green functions for the 
three- dimensional problems is that they must be calculated 


numerically even for the simplest cases if qpiantitative results are 

pn PA "^p 

desired. However some numerical tables are available ^ ^ and 


recently a simplified numerical method based on the Fourier series 
expansion has been given by Mahanty ♦ These Green functions also 
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occur in problems of scattering of band electrons anc^ spin waves 
from substitutional impurities. 

Using equation (2,2) for U|^(k , R ) in .equation (2.9 
the Gieen function matrix elements can be written as, 

,H> ) = - 1.. y (i ia-te' -s ) / , . 

« y, (k )_ tyS > 

k 

, A _ ~ . 

e matrix G is cyclic like the matrix A and hence its elements 

depend upon the difference R - R' . 

Now the equation (2.7 ) can be solved for U . The first 
Bom approximation solution consists of writing Ug for U on the 
light hand, side of equation (2.7 ). Thus the first Bom approximation 
solution of equation (2.7 ) reads , 

2 = Eo - S P JJq ( 2.11 ) 

This shows that U is equal to the incident wave U^ plus 
another part which we wish to represent an outgoing wave for our 
scattering problem. This can be achieved by writing equation (2.10 ) 
in the form , 




with £ ->-0 . 

In the usual way the second Bom apprcximation solution 

will be , 

2 = Ss - S ? 2^ - ( s ^ 2^ 


( 2 . 13 ) 
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am so on. Thus one gets the perturbation series solution. But this 
^oes not show resonance-in the scattering process. Also for lai^e 
portuibations such a series has comrergence problems. 

Instead of equation (2.13 ) we wish to solve equation (2.7 ) 
exactly. In a formal way, the exact solution reads. 


TT rr A A, A A A 1 

U^Uq - GP(i-(-g P)"^ 


u. 


(2.14 ) 


The orthonormality of solutions of the form oi »quHlon(2 .14) has 
been discussed in detail by Klein and others^ We now pir>ceed 
to describe the operational procedure for obtaining an exact 
solution of equation (2.7) like the one given by equation (2.14), 

§ Solution for the scattering amplituda 

In order to obtain an explicit solution we partition TJ in 
the form, 

- “ (scdJ 

where u is a n-dimensionaa vector in the space corresponding 
to the nonvanishing partition ^ of f described in equation 

(2.6). In a similar fashion, the matrix G and the vector V can 
be partitioned into , 


and 


K A 

^ ®( 12 ) ^ 
^(21) ^(22y 




( 2.16) 
( 2.17) 


where g is n x n matrix in P-apaee eorrespondir^ to the 
partition ^ and jio is in space corresponding to a . 
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Let us define the N - dimeneional vectors g, and Sq 

-A 

S = P U 


and So = P a 




(2.18) 

(2.19) 


By multiplying the partitioned fom of P given by equation (2.6 ) 
and using equations (2.15) and (2.17) for the partition of U and 
Uo one gets. 


and 



( 2.20 ) 
( 2.21 ) 


where =* ? u and ^ 



( 2 . 22 ) 


A 

Now multiplying equation (2.7 ) by P on both the sides 


by the left, one gets 


A 

P u 


A 

p u, 




A A A 

- P G P U 


or 


A A 

^ - p G S 

— o ~ 


( 2.23 ) 


which in the partitioned form reads 



whence the reduced equation for ^ 

a “ So - & g s 

so that It = (I + pg )“’^ 
or “^ ( 1 + P g P Ho 



is , 

( 2.24 ) 

( 2.25 ) 
( 2.26 ) 
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Thus the matrix irnrolvef? in equation (2.26) is only of the 

. ■ f\ 

damenflions n x n. For a given case one knows ? and g explicitly. 

But still to find the inverse of(I + ?g)ina stre.ight-forward 

manner is often veiy tedious. To ease this task, one can employ a 

unitary matrix Y vhich block diagonolizes I + p g. One could eYen 

think in terms of a unitojy matrix which completely diagonoli 2 es it, 

but such a matrix is difficult to obtain. On the other hand there 

A 

exists a unitary matrix Y vAxich can be foxond easily by group 
theoretical procedure and which can block diagonolize p and g and 
hence ( I + P g ) into blocks of quite manageable dimensions. Then 
it is very easy to find the inverse of ( I + p g )• The number of 

such blocks, the dimensionality of each block and the number of 

A 

identical blocks for "the block diagonclization by such a Y matrix can 
be known by applying the following theorem from group theory to the 
decomposition of H given by equation (2.4 ) . The theorem^® states 
that if ^ and if p- hag the dimensionality 

V 

then in the block diagonal form there will be n^; identical 
,a X a^ blocks . 


Farther if one writes the block diagonal matrix 



these blocks which corresponds to a certain irreducible representa- 
tion is non-zero, then using each of them one can evaluate 
vectors s^ whose sum will be s . This deconjJosition of s 
is equivalent to projecting the different symmetry adapted functions 
out of s , which belcarg to different irreducible lepiesentAtioruof . 
the point group. This in turn leads to the partial wave expansion 
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of the scattering amplitude. Thus the matrix V found ty g^up 
theoretic procedure not only facilitates the eTaluation of the 
inverse matrix in equation (2.26 ) but also leads to partial wave 
expansion of the scattering amplitude as shown below . The evaluation 
of such matrices is described in the Appendix HI. 

Thus let, 



where T is the blop.k diagonal matrix with h blocks. It io veiy' 
easy to find the inverse of T since this has blocks of very small 
dimensions. Let T be written in the block fom. 



( 2.28 ) 


Now each of the blocks t]|^ , tg, ••• correspond to a unique 
irreducible representation of the point group of the lattice. We 

define the matrices such that has non-vanishing elements 

in the ^ th block only, i^e 

A ' 

1 ( 2.29 ) 

and !rl « 22 (2.30 ) 

Thus (l+pg) «vr”^v 
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^ A 


where ■= V p u, 


y 


( 2.31 ) 


One chtains "fche sgihb rcsul't xf one hloclc ^.xe^onolizes p £s 
A 

well • by the V matrix and writes it ax the sum of matrices like 


/N. 

Mo. ^ LiLiltiplies the ^ th blocks of these anrt defines Sy as 

a,, = T (M^ . 

Th® ^ defined like this belong ti^ particular rows of 
particular irreducible representations of the point greup of the 
lattice-. This identification can be done by identifying the columns 

A 

of 7 - matrix which is done during its evaluation. One would get 
the same result if one finds s by direct inversion in equation 
(2.26) and then applies the group theoretic projection operators 
to s to find the s ^ that belong to given rows of given 


irreducible representations. 


Now going back to equation (2.20) one finds that S can be 
wiutten as , 

h 

S = T' S , 


where 




( 2,32) 




and then equation ( 2.7 ) reads 

h 


2 = a.- ^ 


0 o 


>>=1 




( 2.33 ) 


This equation is the analogue of the familiar partial wave 
analysis result of the theory of potontial scattering. Here the 
number of terms is finite end equals the number of submatrices in 
the block diagonal form of p te-hich in tiirn depends upon the range 



of inteiHction due to the impurity and the point group of the lattice. 

In component foim, equation (2.33) becomes, 

h 

g) =D (k,E) G(a>^. (k^E ) (a.34) 

idieie refers to the lattice sites affected ty the perturbation including 
the impuritj" site and the summation is over all of them. These are 
quite small in number and the non-vanishing elements of refer to these 

only. Using equation (2.10) for G( R,R^) one gets for large R, as 

shown in the Appendix II, 

*1^— ^"bo 

— o 

and consequently equation (2.34) gives 

n(k,E) = U^(k,R) +2Z,g (!£o) SS-fi !So-E)7" (2.35) 

all fe fB|M^ ^ 

where Fo ) is a function of the components of k and k . Here k„ 

is such that a}^(k^) = cD with (k) | having the direction 

of R, For symmetry directions in the lattice ^ represents the 
direction of the scattered wave and often there are not more than one 
JSq satisfying (k^) = ^ along these directions. 

Equation (2.35) gives us the desired solution with the outgoing 

spherical waves characterized by eap(i ^.R)/|,R J \daose co-efficient 

is the scattering anplitude. This point has been discussed in detail 
9 

by Callaway • Thus the scattering amplitude is 

h 
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an(? the desired partial wave expansion can be written as, 

f = 21 V (2-37) ■ 

where f^(k,kQ) = ^ ^ with sumined over n^^ values 

where n^^ is the dinBnsionality (degeneracy) of the th irreducible 
representation . 

For incident waves almg the synimetry axes in the lattice, it is 

quite easy to evaluate the forward scattering amplitudes. Thus one sets 

k “ ISq in equation (2.36) to get the forward scattering amplitude. 

In describing the details of the partial wave analysis ve shall 

follow the classification of S,P,D,F etc. "type of waves according to the 

Wigner's scheme described by Callaway and given in detail in Appendix III, 

Thus having obtained the scattering amplitude one can evaluate the 

differential scattering cross section in the usual way. Also the total 

scattering cross section can be found by integrating over the angles. 

Such a procedure is cumbersome in general, but in the limit of long 

wares it can be done for sinrole cases. In this limit the constant 

frequency surfaces become spherical for the model under discussion and 

hence V cD (k) I lies along k . Then k represents the direction 
k “O ~o . 

of propagation as well as the direction of observation for all directions 
of incident waves which are characterized by k and the angle between k 
and ^ is the scattering angle. 

As mentioned above, one can evaluate the forward scattering 
amplitude for incidence along synmetry directions . For such directions 
of incidence one can evaluate the total scattering cross section 
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by asinj the optical tneorem* That the optical theoiein holds good for 


thcj type of scattering under considera.tion has been discussed by Thoma 
13 15 

and Ludv0.g and Klein • Me feel that the optical theorem given by 
13 

Thoma and Ludwig xdiich can be derived using the results for Q — ( B ) 

9 ^ 

given bc^ Callaway ^ is sufficient for the present work and it is given 

by 


CT' - 

T 


0)^1 

I g^d, , 


(2.38) 

whei'e f(k,k) is the forward scattering amplitude. Here CO^ | . . (k) 
means the second derivatives of cb^i'k) in k space with respect to a ' 
coordinate system that is perpendicular to gradj^c5^(k) and diagonalized 
the tensor j . (k)^ (i, j “ 1,2) and is given ty sgn Ct)^l. (k). 
-dowcver, it must be remarked here that from the experimental point of 


view one is interested not in the total scatteling cross section' ttit'the 
momentum transfer cross section which is required for the calculation of 
themal conductivity. This is discussed in the following section. The 


lind.t of long wa^res is veiy relevant for such calculations because it 


is a plausible approximation for explaining the results for low temperature 
experiments. Also the optical theorem in this limit takes the usual 

k)^for the model under discussion. 

lb is clear that resonant scatterings are included in the above 
calculations. Equation (2.26) shows that _s has det j I + P^j in the 
denominator. Further from equation (2.27) it follows that det | I + p g j 
= det I T-j which in turn is the product of the deteradnants of the blocks 


form of 


417 

nr 


Im 
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' . I A A , 

of ■'’7~ which correspond to different symmetiy types. Thus det | I + P g| 

= det I nT”} = It. b. where b is the deteiminant of the block which 

* J J ;5 

coxrespnnds to a given syinmetiy. These are exhibited explicitly in the 

rexb two chapters. These sub-determinants b. occur as denominators in the 

0 

expressions for f. (k,^) or equivalently one can speak of T. b. 

t \ ' t] ti 

occuring in the denominator of f (k,k ) ,. In general b. are complex within 

u 

the. band, and can be written asb.=b. +ib , where b. and b. . 

3 3,r 2,T 

are read.. Moreover these are functions of the incident energy E. Outside 

the band b .S'O and b. are completely real. Localized modes appear when 
3,1 3 

any of the sub-determinants b^ vanish outside the bend . This can be seen by 
multiplying equation (2.5) by G from the left to get U ^ S* Here 

we do nob have a because we consider outside the band. Now following 
the partitioning described by equations (2.15), (2.16) and (2.6), we get 


^ A 

' g P li 


(2.39) 


This is an eigenvalue, eigenvector problem in p -space, 
equation exists for CO^ determined by the condition det 


Solution of this 

"TTb 


"= 0. These are the localized modes of vibration whose amplitude falis off 


rapidly with distance from the inpurity and these have been discussed in 
detail by many workers^^^^. Since b are Identified into different symmetiy 

J 

types, one speaks of localized modes of different ^yimnetries and degeneracies 


corresponding to the b that vanishes to give rise to these modes. Since 

we evaluate the factors b. explicitly in the present formalism, it becomss 

3_ 

veiy easy to investigate the existence of localized modes. 

Inside the band, the real part b, may vanish for certain values 

Ur 
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of incident eneigy E. It will correspond, to a resonance in the scattering 

amplitude provided that the imaginary part b. . and '^j»r have the 

flE 

same sign and their ratio is small ► This additional requirement comes from 

the consideration of the density of phonon states of the perturbed lattice 

as shown below. We do not use the density of phonon states elsewhere 

althougn- it can be evaluated easily using the above formalism. The 

density of phonon states for the perturbed lattice can be written as 

% ™ ) where the perturbed Green function G» is given bv 

11 


A 

G' 


A A A 

G - G P 


A A A _1 

(I + G P) 


whence the contribution of the impurity to the density of states can be 
written as 


= - i 


('3et (i *GFj)l 

A 


(?..40) 


But due to the partition of P into the only non-vanishing small block 


P; one has 


det j( I + G P )| = det I I + g p I = Ty b^ 


A A 


where we also use the fast that (I + p J ) - (l t J p ) and that the 

transposition does not alter the value of the determinant. Then equation 
(2a40) becomes 

QD. 

J.,, 

■ dE 




... lZ I. f -^1 

T 3 L h J 




db 


3 

3 >r 


_ b 


db. . 


dE 


dE 


^ (b. f + (b f 


(2.41) 
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But at resonance, b =0 and = "I (2,42) 

^ TtL 3,1 J 

is the contribution from a particular impurity state, which w'lien 


compared with the resonance formula 


n (E) = I 


(E-Eq)^ + 


u) 


which at resonance reads 


n(E) 


TT n. 


Ld 


gives the result for the level width 


11 

CO 


b. 

3,i 


bl 


(2.43) 


Here prime denotes the derivative with respeot to E. For resonance 
Ho should be positive. The sme.ller the the sharper is the 

resonance . 


^2.5 Relation idLth. the transport phenomena: 


So far only one defect at a particular lattice site has been 
considered. In actual practice one has a number of such defects randomly 
distributed over the lattice sites. These can even migrate and interpret 
with each other. But at low temperatures and small concentrations, it is 
perhaps not a bad assumption that very little migration of defects takes 
place. We want to consider the effect of scattering under these conditions 
upon the heat transport phenomenan^ but since we will not be interested 
in the eaplicit evaluation of thermal conductivity in the present work, 
the details will be omitted. 

The problem of thermal conductivity can be tackled in two ways. 
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The co-iirentional method is that of writing down the Bolt2anann tran^ort 
.. 36,37 

equation and solving it in the Kinetic theoiy approximation i The 

other raore recent method is through the so called Kubo^® formula. This 

39 

has been used recently by Behera and Deo among others to calculate the 
tnezmal conductivity. However we confine ourselves to a discussion of 
the conventional method. 

In the conventional method one has two distinct approaches which 
give tho same re spit. One is the classical method in which one writes 
the energy flow equation along the lines indicated by Choquard"^^, evaluating 
tho various quantities for the perturbed system. The other is the 
quantum approach through transition probabilities and scattering •ross 
sections. Since in the present formalism one can evaluate tlie scattering 
cross sections, we discuss here this latter approach. One starts wiih a 
Bou.t2mann equation for the rromber of phonons with momentum k and 
studies its time evolution through the equation of the form 

^ ? t y drift ^ 


/coll. 


( 2 . 44 ) 


During thermal conduction, the left hand side is due to the temperature 
gradient and has the form 

dN 

^.grad T 

whereas the right hand side is due to the scattering for the present 
case. In the relaxation time approximation one writes 



( 2 . 45 ) 
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which shows that in absence of the temperature gradient, the deviation 
from equilibrium in the mode k damps out exponentially in the chara ©t eristic 
tine (k) . With this one can solve the Boltzmann equation (2.44) for 
and calculate the thermal conductivity. In terms of heat capacity 
C(k) and the group velocity 3;, tho expression for thermal conductivity is 


K 



C(k) v^ cos^O 


(2.46) 


where 7 is the volume, Q is the angle between y and grad T and the 
integration is ever the first Biillouin zone . Defining the mean free path 


A. (k) = 'X (k) v^one has 


K 


7 

(2 TT )^ 


d^k 


(2.47) 


and y\. (k) is given by^^ 

d 

Aik) 

where is the concentration (small) of the defects, is the 

differential cross sectianat polar angle 0. The weightage (l-cos 0) 
measures the relative change in the con5)onent of the momenctum along the 
initial direction of motion. The quantity 

•= 2Tr 
0 

is sometimes called the momentum transfer cross section. If is 

independent of 0 as is the case for scattering from isotopic impurities, 
when the scattering is entirely S wave type, then equals the total 





(l-cos 0 ) ^)3(©) sin 0 d 9 


(2.49) 


JTTN^ \ (1 COS e) (TpCe) sin 0 d 0 (2*48) 
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cr'.ss section which can also be eTOluated by using optical theorem. Of 
special importance is the Zong wave behavior of y\ (k) . Since tT'uC©) 
cenbains resonance denominators^ the quantity _/V(k) registers a small 
vaZue at the resonance in This in turn implies a smaller value of 

the thecral conductivity K than when resonance is absent. Thus the . 
res'^nance in 0 ^(s) implies a 'dip' in the thermal conductivity curve. 

in interesting feature of equation (2.48) is the interference 
beu'ieen the terms of different syaometries. In particular terms of different 
inversion symmetries, (or parities) namely the ’gerade' and the 'urgerade* 
ternis show interference. This is seen very easily in the limit of long 
waves where one has for example f(k,^) = f + f cos 9 with f_ and f 
independent of 9. Then CJ'ji = f* (ICjICq) f (k,^) 

whence by equation (2.49) 


<rn = ^Tri|fg| 


1 fu 73 - " Vg) 


The second term on the right hand side shows the above mentioned interference. 


This is shown more explicitly in the next chapter- It demonstrates the 

importance of the weightage factor (l-cos 9) in equation (2.48). This 

IT 

point has also been discussed by Erumhansl and Matthew . 

Equation (2.48) shows another important result that for small 


concentrations for which the multiple scatterings from the defects can be 


neg-.ected, the result for defects is just Nj^ times the result for a 
single defect. If one starts from the consideration of TEndom distribution 
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of defects at smll concentrations and averages out the various 
qaantitiesi one aTrires at" this ’"result as has been discassed Klein 
and Vashishta^®. Thus the single impurity problem which appears too 
academic to begin with, is relevant to the physical situation under 
certain circumstances. 



CHiiPTER III 


SIMPLE CUBIC LATTICE 


^ 3,1, The Model ; 

In this chapter we consider a specific model and discuss the 

specific results in an explicit form. This model of a monoatomic simple 

cubic lattice- with equal central and noncentral harmonic force constants 

4 

has been discussed by Montroll and Potts , There are W atoms of mass M 
each. The harmonic force constant Y characterizes both the central and 
non-central forces between the nearest neighbours only. This equality of 
central and non-central .forces enables one to discuss the problem withait 
bringihg in the question of polarization. With lattice constant a and 
Cartesian unit vectors e^, the lattice sites have co-ordinates 


R 


t a % 


+ m a 


-2 


n a 


^3 


(3.1) 


where , m and n are integers. We shall use either R or ( C "to 

specify a lattice site. The cartesian components of displacement from 

equilibrium of the atom at site (L >m,n) are given by X i Vj and 

“n ' t mn 
4 

^ tmn* discussed by Mont roll and Potts , one can think of any of 

these X/ for example as a parameter associated with each lattice point 
tmn 

( (,,m,n) and then this model is equivalent to a scaler model. With 
this it is possible to write the e(jiations of motion of the perfect 
lattice in the form of equation (2»l) 




0 


(3.1) 


/\ 

where A has the matrix elements given by 


, 28 .. 



* «29 » • 


A(lmn, ^ is,m' ^ n,n' 

^ m,m’ ^ n,n'^^l+l, t! ^ ^^m+l,in' 

^ rsr-lf ^ ^n+l,n' ^n-1, ^ (3«2) 


and 


mn 


y (k, (,mn) “ 

The same matrix element A((tmn, t^m’n’) will describe the motion when 

U (k, ( mn) are chosen to be Y/ or Zt 
~ Lmn '^mn 

The plane mvB solutions cone ponding to equation (2.2) are 
(k, l^mn) = exp jj. (k^L + k^m + IJgn) J (3.3) 


which correspond to the dispersion law 


n O ^ 

CO (s = — (5 - 


(3,4) 


Now let the defect atom be situated at the origin of the ^foordinate 
system. It has a mass Mt- M and is ccsupled to its nearest neighbours ’ 

A 

through force constant V + aV- The pertuibation matrix P has 
non-vanishing elements only when the indices refer to the impurity site and 
its nearest neighbours. D<^noting the nearest neighbour vectors by 

A 

the non-vanishing elements of P are 

P (o^s) ■= am J - 6 4 V , P (^, ^ V 


and 


P (Qir^) “ P (Pft.O} = ^ y 


(3.5) 
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Here ^ can have eompocents (+1> 0,0), (0,+!, O) and (0,0,+l). 

Q 3 .2 . The partial wave .analysisi 

^ Is descrihed in Appendix III, the reducible representation of the 

cubic group 0^ in the space of the affected lattice sites ( p -space) 

42 

reduces in the usual notation to 


p. 


2A, + P + E 

ig lu g 


(3.6) 


A 


This leads to the unitaiy matrix V given explicitly in Appendix . 

s 

III. Now using the procedure outlined in the last chapter one can find 
the scattering anplitude. Here we follow the notation for the Green 
function integrals described in Appendix I for the simple cubic lattice. 
Also the matrices^ and ^ are explicitly written with same indices for rows 
and columns as ere used to label the rows and columns of in Appendix HI. 

With these the elements of block diagOT*! form of (l+P^ > namely the 


t matrices of equation (2 #28) are given by 




t = 
1 


s 


A A 

^3 


% = ^6 



-6‘^H 


- 6^(c„- E^) 


lA , +AC^ + as 
Ai 0 ° 


l/E. 


^■D 


(3.7) 


where 


A 

H 

Q 


IE, _ 6 

AY 


AE 


+ ^41(1,1, 0,Eg) + I(2,0,0,Sg)^ 


o o 2Y ^ 

1 + 6E_ -- C„ + ^ 4l(l,l,0,Eg) + l(2,0,0,Eg)'| (3*8) 


AV 


Vl 


J 



. . 31 ,, 


with E 


^ I (1,0,0,E ) 




and 


and 


i 


q[i + Ai{Uc^ - 

|- e4Y h[(VV] 

(3.9) 

1 

)I (0,0,0,E ) - 

■ I (2,0,0,E )7 

(3.10) 

;iY 

L 

J 



> I (0,0,0,Ej + I (2,0,0,E„) -21(1,1,0, Ejl 

z-/ ' 

L 


J 


Using these one can evaluate the scattering amplitude as discussed in 
the last chapter and the result is that the scattering anplitude 
f(k,l^) can be written as the sum of three terms 


where the S-like wave amplitude is 




go(l^o) 


aMcO^ -6 AY + ? 4 l(l,l,0,Eg) 

2 Y C 


+ 1(2,0, 0,E^) + I(0,0,0,E_) 




1> ,.Al?iL I<1,0,0, 


E, 


X^^(coa k. . 00 s I(0,0,0,E^)J^ 

X^lgooak^^g 

^oi 3 ] (3.13) 


The P~like wave amplitude is 


.!S|^(AV) g sin k^ sin ^ 0 ;^ (3 J-4) 



. .32 


and the D-like wave amplitude is 






3 F. 


D 


3(cos kg. - cos k2)(cos k^^-cos k^^) 


+ (2cos kj^” cos kg “ cos cos k cos k^g - cos k^^) 


(3.15) 


When the substitutional impurity is an isotope, the ? and D- wave 
amplitudes vanish and the S-wave amplitude remains in the foim 




gs(ko' 


F 

Si 


where F 


Si 


1 - ^ 1(0, 0,0, Eg) 


2 'i 


(3.16) 


(3.17) 


8 


This expression has been obtained by Callaway. 


3.3. The scattering cross section, its long wavelength limit 
and resonances : 


Here k specifies the incident direction. Is already discussed, 
along the ^jnnmetiy axes of the ciystal, coincides with the direction of 
scattering due to the behaviour of the constant frequencgr surface along 
these directions. This facilitates the use of optical theorem to evelua'fce 
the total cross section when the phonon is incident along cne of the 
^mmetiy axes. JO.so in the limit of long waves, this can be done for 
any incident direction. The scattering angle in this case is that 

between k and k . ' 

— —0 

To demonstrate the use of optical theorem, we evaluate the total 
scattering cross section <5^ for phonons incident along (1,0, 0) axis 
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using equation (2.38). Then k = (k,0,0) and we set ^ = (k^,0,0) to 
evaluate the forward scattering amplitude, k and k are equal in 

■" —Q 

magnitude. With k = (k,0,0), g (ko) = 2^®”^ scattering cross 

o S ^ ^ J 

2 

section can be expressed in units of a where a is the lattice constant^ 


^ ^ ^ - ^ ( Im D ReFs - Re D Im Pg) 

" V-nk[|Fs|2 


^42- (Im Fp) sin^ k f - 4V (cos k-l) (Im F ) j (3.18) 

12 ^ . 1 ^ l2 ^ i 


where 


D = AM 60^ -6 A.V+ 4 M of 5^4 1(1, 1^0, E.) + l(2,0,0,E^) 

2 V L 

+ I(0,0,0,E^)J + 4 ^1 - .. 4^ - ^ . - l(l,0,0,Eg) j. 

_ -I ^1 - I(0,0,0,E3)J (2+ cos k)! (2+ cos k) (3.19) 


l(0,0,0,E ) y (2+ cos k) (2+ cos k) (3,19) 


^ ^ ^ S(0,0,0,Eg) + ^73(0,0,0, Eg) 

+ 12 C(1,0,0,E ) - 4 S(1,1,0,E ) - S(2,0,0,E^)?+ Ak.(3-E ) 

2^ M 

X , je ^C^(l,0,0,Eg) S^(l,0,0,Eg^ C(0,0,0,Eg)^4C(l,l,0,Eg) 

+ C(2,0,0,Eg) - C( 0,0,0, Eg )J + S(0,0,0,E^) ^4S(l,l,0,Sg) 


+ S(2,0,0,E^) - S(0,0,0,E^) 


}] 


(3.20) 



.. 34 .. 


Im (3-E„) C(0,0,0,E„) + 7C(0,0,0,Eg) - 12S(l,0,0,Es) 

M ® 2 y 



- 4C(l,l,0,Eg) - C (2,0,0,E„)J. + 


A 


M 


2V 


M 


(3-EJ 


He 


In Fn 


X 1^-120(1,0,0, Eg) S(1,0,0,E3) + C(0,0,0,Es) ^4S(l,l,0,Eg) 

+ S(2,0,0,Eg) - S(0,0,0,Eg)J + S(0,0,0,Eg 
+ C(2,0,0,E3) - C(0,0,0,Eg) j J (3.21) 

^ ^ ^ 2^=^ S(o,0,0,E3) + S(2,0,0,E3)^ (3.22) 

= _4i 


Re Fp = 


In Fp = 


2 y ^C(o,0,0,E3) + C(2,0,0,E3) J (3.23) 

1 + -||_^S(o, 0,0,E3) ^ 23(1,1,0, Eg) - 3(2,0, 0,B3)j (3.24) 
~i-|^C(0,0,0,Es) + 20(1,1, 0,Bg) - 0(2, 0,0, Eg) ^ (3.25) 


Here we have used the 0(p,q,r,Eg) and 3(p,,q,r,E ) defined in Appendix I 

28,32 ■ 

and their values have been tabulated. 


In the special case when the substitutional impuri'ty is an 
isotope , 0 so that D = A and F . = 1 - /^ ^ l(0 

^ 2 V y 

which leads to the result 

c57i ' 7 


,0,0,Eg) 


1 / Am \ 

sin l4 2V J 


o) 


4 

We define 0)' 


oS 


uL 

M 


^3i 


2 - Iin(l (0,0,0, Eg)) (3.26) 

(3.27) 


and in the limLt of long waves 

and then with sin kHkk, equation (3.26) yields 


47rMa]^£ 


(3,28) 
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<r, 


1 _ _ f Am >,2 £;i>^ 


TT 


M 


^!s '"sil 


in units of a, (3.29) 


This is exactly the expression obtained ty Callaway . The scattering is 
purely S-wave type . 

The long wavelength limit of CTJ of equation (3.18) is given by 


' Y"' [ I Fsp 


Im D Re Fg - Re D Im Fs 


] 


2 k^(lm F^) + k^(-lBi ^) 


31 Fp 


(3,30) 


where the quantity D now becomes 


D - AMcO^ -(k^ A) 41(1,1, 0,Eg) + I(2,0,0,S 

2 V ^ 


+ 71(0,0,0, Eg) - 121(1, 0,0,Eg)J+ ^I(1,0,0,E 


y 

- I(0,0,0,EJ 


} 


(3.31) 


where l(p,q,r,Eg) are to be evaluated in the limit of long waves. 

4 

Retaining terms upto the order of W ^ the D--wave amplitude does not 

6 

contribute because it is of the order of c5 and hence 

n- . _1_ Yi , b _4V f 


fAlY _i_l (AL Y 

^ Y / hpl^J V <^os' 


(3.32) 


in units of a „ where the constant b This result shows that 

? 


in 


the long wavelength (low frequency ) limit, the scattering cross secticai 
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is proportional to cO showing Rayleigh type of scatteilng, provided 

that there is no resonance which is seen ty examining Re Fg and Re Fp . 

The S-wRve part of the cross section is affected ty both A. M and AV 

whereas the ?-wave part depends only upon A"]/* 

It is worth while to point out that the number of partial waves 

contributing to the total scattering cross-section depends strongly upon 

the direction of incidence. This is discussed in some detail in the 

next chapter. For example if we evaluate the scattering cro.-vs section 

for the ( 1 , 1 , 1 ) direction of incidence the D-wave part will not contribute 

to an expression like that given in equation (3.18). 

A resonance in the scattering cross section will occur for that 

value of the frequency at which any of the quantities ReF 3 , Re Fp or 

Re. Fj) vanishes and the width = EAZ is small and positive. 

^ (Re F)’ 

The real and imaginary parts of Fg, Fp and Fj, are given by eqviations 
(3.20) to (3.25). One can study the resonances by taking given values 
of - ^ M and . ^ and examining the criteria for resonance. The 

M y 

figures 1 and 2 give the results of such a study. In these the value 
of the parameter E , defined in Appendix I is found at resonance for given 


values of 


A M 


and ^ V . ¥e find that the ? and D resonances depend 
M Y 


upon 47 only. For sufficient decrease in force constant the P and D 
resonances appear at the low frxsquency side of the band.- If this is 
increased the resonances appear towards the high frequency side of the 
band. Indeed in figure 2 if one changes the sign of 

corresponding Egj^ simply reverses the sign having its magnitude intact. 

The S-wave resonances depend upon Am as well as A V • ® 

large increase in mass, resonances appear towards the low frequency 





.. 39 . . 


(large positive S^) side of the hand, no matter what the change in 

force constant is. On the other hand, for small changes in mass, the 

change in force constant plays a dominant role in determining. the resonances. 

4 

At resonances, the scattering cross section does not have a ---gO 
dependence in the limit of long waves. Since- the real, parts" of Pg.and 
Fp become zero and the imaginaiy parts of Fg and Fp are. proportional t.o 
C^J , then equation (3.30) for ■ 0 ^ shows ' that -it is proportional 
"fco s’t resonance.- Further the xd-dth- of the ..resonance is 

proportional J * Thus the value of at TOsonaRce beoonBs- 

leige .by a fa-ctor xri-th- small Ci) • Using the expressions given 

by equations (3.20)., (3.22) and (3.24), oJ*fe can study the localized modes 
of -S ,P and F symme.trie s . 

3-,4 Phonon scattering and the mean free path ; 

We evaluate the mean free path of phonons in -the limit of long 

^ '4* • 

waves hyusihg equation (2-^48). The differential sjiettering -c»oas .section 
Q this limit is simply tho squa-re of the scattering, an^^tudo 

which is given by 






0 ) 


0 ) 


p 

(3.33) 


oS 


with b fid 1.14 and © is the scattering angle. 

Then 

T 

^ = 2 IT Tj f(©) I (l-cos e) sin © d 6 




0 


♦ • 4:0 ♦ • 



Tho last 'berm on tho lO-ght hand side is the interesting toim showing 
tho interforonco botwoon tonns of opposits 'paritios' . Here 'parity' 
signifies the behavior under inversion at the centre of inversion. The 
S-type of teim does net change sign under inversion and it is the so 
called 'gerade' term and the ?-type is'ungerade*. 

At resonance the cd ^ proportionality of J\_ 
hold and the dominant term can be shown to be proportional to Cd 
as in the case of 

Thus the simple model considered here shows the important 
features of the scattering of phonons except the question of polarization 
which is discussed in chapter V, The long wave length Rayleigh type 
of scattering end the modifications at resonances are explicitly 
obtained here. The directional dependence of the scattering cross 
section, the interference of terms of different ^jrmme tries in the mean fie e 
path expression are the other in®)ortant features of the J^rocess and 
these are further elaborated in chapters IV and V. In the following 
section we consider a model of diatomic simple cubic lattice to 
demonstrate the ease vdth which the results for the mono-atomic lattices 
Can be carried over to the study of diatomic lattice* In fact this feature 
adds to the useftilness of the results for mono-atomic lattices. 


Ca) does not 
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3,5 Diatomic Simple cubic lattice 

Let a diatomic simple cubic lattice be considered idiose 
atoms of type I have mass each and those of tjrpe II have mass 
each. The central and non-central harmonic force constant 
between nearest neighbours is denoted by Y . One can write down 
the equations of motion and perform the so called M transformation 
used by Montroll and Potts^“^ so that the equations of motion for 
the diatomic lattice becomes identical to those of monoatomic 
lattice vath M C() in place of M viiere, 

M* •= 6 Y + ^ si) - & i) (3.35 ) 

Then the dispersion becomes 

M* £*)^(k) = ey - 2y 2 oos. (3.36 ) 


and this leads to 

\/ (Mj^ cO ^(k ) - 6 y ) (M2£d^(k ) - 6 i ) ^ Z ^ ^ cos 

(3.37 ) 

Solving this equation for ^(k ) one gets the optical 
and acoustic .branches as. 




^ ) r 1+1 


and 


4 ( 1 - i cos k^ ) 1 

)" _ 

(3.38a) 


u)^ /i- 

ac ^ 


4 (l - + cos k, ) 


M^M3_(. 




J 


(3.38b) 


The maximum frequency in the optical band is d)2 


ipM 


= 6 •y (-^ 1 ^ + ^ ndnimum frequency in the acoustic band is 
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cO ^ = 0 . The mininiam in the optical band Cx) and the 

a cm opm 

TPaxiinum in the acoustic band Cl)' ^ lecend upon the relative 


magnitude of and . For > Mg and 

whereas for M^^ and c5 ^ = -n--- • The gap 

1» opm Mi acM 


6 l > , ^^2 6 7 

^2 acM 3^ 


between the two bands is of 


magnitude 6 ( 1/M^- I 


In analogy with the monoatomic case the Green function integral 


becomes , 


G ( m , n ) = - l/N 


exp ^ i(lcj, i + k^m + k^n ) J 


(i(^t 


2 X 1 L JJ 

"s y“N ^ T”'l; cos'k" 


+ kgffi + k^n 


(3.39 ) 


where E = 


(M^ d)^ - )(M^a)^ - 6 y ) 


^0)2 

acM ^ 


cO 2 


(3.40 ) 


(3.41 ) 


It may be noted that at each step one can check. the results with 


that for the monoatomic lattice by setting ^ “ ^2 ' quantity E 

6^2 

goes over to E = 3 - for M, = M„ . In terms of the E defined 

0)2 12 
OS 

above , 3 E ^ 0 defines the acoustic band and 0^ E ^ - 3 

defines the optical band. Within the gap, B has an imaginary value. 
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Thus taking appropriate values of E one can consider the 
scattering of phonons in both the acoustic and optical bands. Optical 
phonons do not contribute to the low temperature thermal conductivityj 
but it is interesting to consider the scattering of optical phonons from 
the substitutional impurity. 

Consider a substitutional impurity put at a site of an atom of 
type I x’hich may be taken to be the origin of the coordinate system. 

Let its mass be + A M and the force constant be "Y + 

The elements of the perturbation matrix after the transfonoation 
get modified into the following, 


P (0,0 ) = ^ ( - 6Y ), P (0 , ) = p (^,0 ) = Ay. 

/I y 


P (R , R ) = 
-n ’ -Ti ^ 




Following the method described in this chapter ^ after performing the 
M transformation) one gets the following resonance denominators also 


obtained by Mitani and Takeno 
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.iAl. d)la,o,o ) 




A M 


(3.42) 


Tp =1 


and Ft) = 1 - 


V 1 

( 0,0,0 ) - G ( 44)^,2,0,0 )^ (3,43 ) 

1^,0, 0,0) + G ( 0)^,2, 0,0) J- 



- 2G(CO ,1,1,0) 


where 5 - oy) / (M^ - 6 y ) 


(3,44) 
(3.45 ) 


'tithin the bands, the parameter E is real and the Green functions 
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occuring above can be written in terms of I ( P,q,r, E ) as, 

G- ( P,q,r,) = - ~~ I ( o,q,r,E ) 

which in turn can be written in terms of the tabulated functions 
G (p,Q,r,E) and S (p,q,r,E) . Here cO * is related to E by the 


(3.46 ) 


Lated to E by the relations. 


(i) Optical band. 


‘ L' o?m acM ^ v/ '' ^opm acM 

+ 4 (eV9) Cd^ ~j 

opm acM — i 

and (ii) Acoustic Bahd, 


CO = 


i ( gO 


opm acM- ^ opm acM 


.2 .2 


( 3.47 ) 


+ 4 (eVs) a>^ 

opm acM 


(3.48 ) 


Using the above expressions for the resonance denominators 
fiitani and Takeno ^ have investigated the localized modes both above 
the top of the optical band and within the gao between the acoustic 
and optical bands. We have investigated the resonances in the acoustic 
and optical bands and the results are given in the figures 3 to 6 for 
the cases; (i) ^ ^ ^ ^ , p = 

= -1- . These are the same as those used by Mitani and Takeno^^ for 
investigating the localized modes. 

The change in mass does not affect the resonances in P and 
D like waves. Resonances only in the optical band are shown for 
these waves. The result is very similar to that of the monoatomic 
lattice case and as the force constant increases, the resona.nce 
frequency moves to-wards the high frequency edge of the band. 
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Figure 5: Typical S-wave resonances in Acoustic Band of diatomic simple cubic lattice * 
Curves a, b, c, and d are for Mj > p=l /2 with -.5, -.3,+*! 

respectively and curves a’, b’, c' and d^ are for M 2 > p=2, with Ay/y 

= -.5, -.3, 0 and +. 1 respectively. 



Figure 6: Typical S~wave Resonances in the Optical band of a diatomic simple cubic lattice, 
lattice. Curves a, b, c and d are for M|^> M 2 , p= 1/2 with ^ 

0 and 4*. 1 respectively and curves a‘, b’, c* and d' are for M^> , p”2, with 

-.3, 0 and 4*. 1 respectively. 
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In all "bhe fignrss E denoles Iho value of E at resonance * For 

ri 

the case which implies that the substitutional impuiity has 

been put at the site of an atom of mass (heavier mass), comparatively 
low values of increase in force constant produce P and D wave resona- 
nces in the optical band. For higher values of increase in the force 
constant one gets these resonances if an atom of mass M is substituted 
by the impurity . 

The S - wave resonances are interesting. Tn the acoustic band 
one gets the usual result that as the increase in mass becomes large, 
one gets low frequency resonances, no matter what the value of the change 
in force constant is. This remains true for both the cases of (i) 
and (ii) . But in the optical band, one gets different behaviors 

for these two cases. I'Jhen the site of the heavier atom is substituted 

Mg) , one gets resonances for a decrease in mass and the smaller 
masses give rise to resonances nearer to the bottom of the optical band. 
'Aien the site of the lighter atom is substituted, one gets resonances 
for which the smaller masses give rise tn resonances near the top of the 
optical band> It is interesting to note that from E=OtoE = -l 
large decrease in force constant is required to get resonances ( not 
shown in the figure). This can perhaps be attributed to the behavior 
of the density of phonon states in this interval. 

From the above analysis it becomes clear that while in the 
monoatomic case, only heavy masses can produces resonances and light 
masses give rise to localized modes, in the diatomic lattice light masses 
can produce localized modes above the optical band and resonances in 
the optical band. In addition one can get localized modes in 1316 gap. 

f 

The present analysis, along with that of Mitani and Takeno'^^ leads to a 
rather complete analysis of the problem for diatomic simple cubic lattice. 




CHiPTEE IV 

BODY CENTEED CUHTC Af'ID FACE CENTEED CUBIC LATTl CES 
^ 4 «1 The Model 

In this chapter the results of a study of the problem of 
scattering of phonons from a substitutional impurity in a simiOlified 
model of face centred cubic and body centred cubic lattices are 
presented. The method of analysis is the same as that presented in 
chapter II. The displacements of the atoms in the lattice are assumed 
to be scalar parameters idiich depend upon lattice sites. This is not 
realistic for lattice d.vngTnics as there would be no polarisation of 
phonons in this model. However, sodb essential features of the process 
of scattering are brought out even in this model and one can get a 
qualitative insight into the problem of scattering in more realistic 
lattices. As shown in chapter V, the results of the present chapter 
can be obtained by a sditable choice of force constants in a realistic 
model. The present analysis is, however, oompletely aoplicable in the 

problem of scattering of spin waves or band electrons from substitutional 
impurities. 

In the two lattice models considered here each atom has 
mass M and nearest neighbour harmonic force constant , except 
for the impurity atom which has mass M + AM and is coupled to its 
nearest neighbmrs through a force constant Y + Y • The impurity 
is taken to be at the origin end. the lattice vectors E have 
Cartesian components'^, 

^ (t + m/2)a ej + ( m/S )a gg + (n + (4.1a) 
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cr;.d 

iCL'^ n)a £]_ + -^C i + m)a e, + ^(m + n)a . ( 4.1b) 

rier=5 e is the length of the side of the unit coll and e£ are the 
unit vectors along the Cartesian axes, and , m, and n are integers. 
The total number of atoms is N. The subscripts 3 and 3 refer to 
B'ody Centred Cubic and face centred cubic lattices respectively . 

The frequencies (squared) are given by, 


CO^(k ) >= (l *“ cos ■ 

. OB ? 


\ ^3 

cos — ^ cos — — ) , (4,2a) 
2 2 ' ^ 


and 


2 

P 


OF 


k 

2 


Co (^)«=!Co (3 — cos ‘ — 1— cos — cos — ~ cos 


2 


k 

/* 

< 

2 


2 


kg k 

_ cos cos — i- ) , ( 4,2b) 

2 2 ' ' 


sxan 


vhere to =8 ”//M, CO = 4 /M, and k. are the Carte 

OB ■’op 1 

components of the \rave number k. These can be obtained ty writing 

the equations of motion in the form of equation (2.1 ) with the 

A 

components of A given by , 

B,F 

^ ^ ^ ^ ® .^m*m ^n^n ~ Sn^n^ ,mtl ^m’,in-4^ 

^ n^n^ (^+1 ^ m* ,m-l (_-! ^m‘'>m+l ^ 

ft * ^ ^ n*n+l C m' ,nH-l £ n* ,r>-l ^ 

^ ^ 

bm* bn*yn+l 

( -1 ^ () 

S7lfl 


(4,3a ) 
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A ( ivarif t’m'nO = yri2 ^i,i C, O ^ -( iC . (S , +)s ) 

^ L om'in dn^n f * m’lii'' ^n’n+l 

) ? D )) A J 

^ C’ f ^n’ n ^ ^m^m+1 ^ m’lQ-l^ ~ ^ m' m ^ * 1+1 ^ ^ h /— 

>> A V A ) } j 

~ ^ t. U ^ ^m' m+l ^ n~l^ ra-l ^ n^n+1 ^m f m^ » (+1 C n-1 


^('(-1 S': 


n*n+l 


^n;„(<r(;/ +1 ^ m*m--l ^ S^m^m +1 ^ 1 

' • -^ ^ j ^ — i 


(4.3 b) 


The pertuiba'fcion inabrix P of equation ( 2.5 ) for "these cases 
has non— Vanishing elements only when the indices refer to the impurity 
site and its nearest neighbours. ¥e shall denote the nearest neighbour 
vectors by . The non-vanishing elements for the two lattices are, 


P (S,2) = AK co^ -(Zg^j) aY, p (S„,4) =-aY > 

p (a.s^) -p (£^,2) = aY , ( 4.4 ) 

where Z = 8 , and Z = 12. 

o F 

The unpe rtuxbed nome.l modes are given by, 

*= IT^ exp (i k.g/e )'■ , ( 4,5 ) 

The elements of "the Green function matrix, along with 
their low frequency expansions are given in Appendix I, 

^ 4.»2 The Partial wave analysis 

Here the point group of the lattice -is the cubic group 
The representation of 0^ generated by the displacements of the 
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affected lattice sites reduces, in the usual notation'^ to, 

= + A , (4.6 a) 

1 B ig lu 2g 2u 


and 

= 2A +F +E+F +F . (4,6 b) 

i F Ig lu g 2g 2a 

One can use the method of projection opeiB.tors to obtain the 

A 

symmetiy vectors which constitute the columns of V matrix as 

A /\ 

described in Appendix III, where V and Y are found e3cplicitly. 

B ^ /\ 

The indices of rows end columns of the 7 matrices as 

A 

specified in Appendix III for 7 and 7 are used similarly to 

A A B F 

A A A 

wilte the matrices Pg, Pj, and G^, in their explicit forms. 

A . . . / 

The inverse block matrices t^ occuring in equation (2.28) 
are as follows. For the body centred cubic lattice. 


t 


IB F 


SB 


-8-^A'/(gQg - gj^g) 1 + Aj 5 gpg + 8 A /g. 

A A A 1 


IB 


AAA 2. 

-5 • t = t = . t t *= t 

2B 3B 4B Fpg " 5B 6B 7B 


T ; 

IB 


and 


Here, 


8B p; 


FB 


A_ *= 4M - 8 
B 

• Ab glB + a/ ( 3 ggB * JSjB * Sob ^ ' 

- (1 am (1* Ap ggj *BAVs^) 

- 8 Bg a/ (ggj - gig) - 


(4.7 ) 

(4.8 ) 
(4.9 ) 


(4ao) 
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PB 


" 1 - A'/(gQg + ggg - §32 - , 


F =1 
DB 


“ ^2B " Sb ”^^46^ ^ 


and^ 


- A -/ U 


- M (gQg - Sggg + 3g^^ - §45 ) • 


(4.11) 

(4.12) 

(4.13) 


The notr?.tion for the Green functions g is (ie scribed in 

iB 

Appendix I. With these^ the scattering amplitude is, 

’ K'> = "'sB fpB ^o) " 


where , 


- 

^b^^^ r 2 

^SB = -it: L (1 + am a) g,, - m„) (A. 


(4.14) 


•SB 


AB 


B 


" S - Sft/mj (1 _ C^C^Cj) + 8 0^,0^jC^3. Af X 

1^(1 - 0^0383) (1 . * 8 a / ) 


OB “IB" '“B ■ ' ”l“2”3" J ^ 


+ S„C_C,S „G _G , + S_G^C^S ,C ,C ^ 

2 3 1 o2 o3 ol 3 1 2 o3 ol o2 


f (k> IL) 

DB “O 




] 

[=^^=01^2^03 " ‘^2^¥o2=o3Sc1 


(4.16) 


+ C_S^S^C -S ,S ^ , . 

3 1 2 o3 ol o2 f * 


1 


and. 




(4.17 ) 
(4.18 ) 
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Here C = cos (k./2), S. = sin (k./2), C . = cos (k /2) anf^. 

^ 1 ^ x' 01 or 

Sq£ = sin (k /2). The expression for g^Ck ) is given in lppenf?ix II, 

OX a 

Similarly, for the face centred cubic lattice one obtains, 


^ 1 Am 4a)\ 


= -1 . 


IF 


12^^V(gop - g^p) 1 + + 12Ayg, 


F'^OF 


A A A 1 A 

^2F " ^3F " ‘‘^4F " ' ■^5F 


'6F F ^ 

or 


A A /\ 2 A A A 2 

^7F " ■^8F " "^PF " ~-r ^ ^inF = "^11 F "^IPF^-V 


lOF “ %1F =^12F"T^" 


where, 


Ap = AMcO^ - 12 A y , 

mp = A^^ „ +AV(g^ + g + 2g + 4g + 4g ) 
V f°1F ^ • cF 5F 2F if 4F 


(4.20) 

(4.21) 


F = (1 + AM ) (1 + + 12 AVg^p ) 


12 ^^OF ~ ®1F ^ 


(4.22 ) 


F _ = 1 + A 7^ - g . - 2g + 2g + g ) , 

FF 0^- IF 4F 5F 


(4.23 ) 


1 gQp + ^^IP *' ^^2F ^ ^®4F ^ %F ^ * (4.24) 




(4.25) 


and , 


^FF ^ ■^’Ay (— gQp + ^^IP “ ^4F ^5F ^ * 


(4.26 ) 


The scattering amplitude in this case can be written in the form, 

fl.(k ,fc„ ) = ) + f,, (k,k_ ) + (k , ) 


■SF '^“0 


?F '-"-o 


+ (k f (k >k ) , 

DF2 “o FF - -o 


(4.27) 
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where^ 




^SF 


(l +AM (X) - Eip) Sip + 4:A~X C^C 




^r2 


+ CgCj + C3C3^)J _ 4Z57'inj,{3 - CjCj - C^Cj - CjCj^) 

+ 4a>f ( ) X 

r VoF ^ (3 - C^C^ - ) 

- (goF - hr> D‘f * ^^VCCjCa + CaCj + CjC^ )j jl (i.ss ) 

Sl(®2 %)Sol(Co2 * <=03^ 




fFp(i gp(4) 


+ 82(03 + C^)Sq 2 (^o 3 ^ ^ol ^ 

" ®3 «=1 * Cs )So3('=ol" “o2> ] 


(4.29 ) 


^DFl(“^^) “3" F^ 


DFl 


«f(4.) [ 3(^=1 - 


+(20203 - 030^ - G^02)( 20 ^ 2 ^ 03 “ ^o 3 ^ol" % l '^ o 2^1 


^DF2^~ F^) 


“ P^“ j^®1^2^ol^o2 ^2^^o2^o3 


+ S-S- s ^ 

3 1 o3 ol 


and; 


j . ( 4.31 ) 

^pp(k ,k^) « - gp (k^)£s^ (O2 - 03)Sq^(C^2- ^^ 3) 


+ 82(03 - 0^ )-o2^^o3 ” 

^ ~ ^2)"o3(^ol - ^o2> ] • 

Here we have distinguished between the two D-wave 


(4.32) 


amplitudes and f’pp2 "tkat arise due to the two imeduoible 

representations E and F. of 0, . This identification is described 
g 2g ■ k 
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in the cheracter table given in Anneni^iy III. It nay be mtef^ 
that above is very similar t^ for the simple cubic 

lattice given by equation (3.15). This nay be expectr-t c^.ue to the 
fact that bnth of them arise fue to the irret^ucible representation 
of . At the same time, is very similar tn both 

arising from the irreducible representation F 

2g 

For the isotopic impurity case only the S-wave amplitude 
remains in both cases, and the form of the amplitude is the same 
as that in the simple cubic lattice given by equations (3.16) and 
(3.17) . 

^ scatt erin g cros s-section, its long w pvelength limit and. 

resonances 

As EBntioned in the last chapter, k specifies the incident 
direction and k^, the direction of scattering viienever k^ is 
along any of the ssnnmetry axes. The total scattering cross-section 
is found by using the Optical TheoTOm. An interesting feature of 
the scattering process is that the number of partial waves that 
contribute to the total scattering cross-section depends on the 
direction of incidence of the phonon. This is expected, in view 
of the fact that the 'potential' due to the impurity is a non- 
central one, and its 'range' depends on the direction in which 
it is viewed. ¥e will discuss the results for the total scattering 
cross-sections corresponding to a few incident directions to 
demonstrate this effect. 

The forward scattering amplitude, xjhich is necessaiy for 
the use of the Optical Theorem to evaluate the total scattering 
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cross secti-r, is epsily found when the- ^irectio'^ of incidence is 
(100), (no) or (ill). One can then set k = . For incidence 

along (lOO) direction, 


k = (k,0,0,) 

and in the Body Centred. Cubic lattice, 

1 


and 


(k/2). 


^(100) ^ _1 f 1 

TB 4Ysrn| 


(Jni Dg^)(lle Fgg ) 


Dg^)(lm sin^(k/2) .(im F^^) 


P?bP 


(4^33) 

in units of a^. Here Dg^ is given ly, 

2 

Bbi = g^g k/2) -a4i'}^g(l - cos k/2) 

+ 8^y cos(k/2) ^(1 - cos k/2 )(1 + Agg^g + S/i^g^g) 

- (SqB " (Ag + B^Vcos k/2 ) . (4.34 ) 

shovjs no D and F wave contributions. For (ill) direction. 


C~iioo) 

TB 


k = k^' = (k/ 32 , k/3^/ k/3^ ), 


and , 

(5^ 


(ill) 


m-n "*""''"7"' — 

27(12)'^ Sin(k/122) Cos^(k/l2'^) 
^ ~ ^B2 


__1 

I'W 


24aV 

PPBi^ 


)(lm Fgg) j 


cos^(k/l22) sin^(k/l22') (im Fpg) 


- cos^(l5/l2'^) sin^(k/l2'^) (im F^g ) 

8in®(kA2’^) (Im Fj^) *j , 


1 ^BB : 


^fbP 


(4.35 ) 
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in units of a .Here is given by, 

Dg2 = (1 + AM - mg)(Ag + c'AV’cos^ k/l2^) 

- - cos^ k/12'^ ) + Six\oos^ k/l2'^ X 

jjl - cosVl2^) (1 + Agggg + SAVg^B ^ 

“ cos'^ k/l2'^ ) X . (4,36 ) 

Thus shows contributions from all the partial waves. 

For (no) direction, it can be shown that the F wave 
contribution is suppressed. 

In the face centred cubic lattice the corresponding 
expressions for the total scattering cross section can be found 
in the same way,. For example. 


(T'Ciio) 

TF 


2V tsMsin k/^)(l + 003 k/8S) 




^(Im Dp) ( Ee Fgp ) - ( Re Dg ) (im F^^ )^ 
sin\/82 (1 + cos k/8~)^(lm F^^) 
cos^k/8‘5' (1 _ cos k/8^)^ (im F ) 


tF |2 .vw ^DPl- 

' DFl’ 

4 A V . 4 

sxn k/8 (im F ) 

F DF2 ^ 

' DF2' 

4 A "Y pi- X o 

^ fZ — is” (1 - cos k/S^) (Im Fpp ) 

* ^FF* 


(4.37) 


in units of a > with given by, 


By = (1 + w^g^p - mp) H- 4AV(2 cos k/8^ + cos\/8’^)J 
— 4AV®f(3 - 2 cos k/8^ - cos^k/8^) + 4aV(2 cos k/8^+cos^k/8^) K 
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(1 + ApgQp + 12^gij,)(3 - 2 cos k/82 - cosVe-) 

- (gQ^ - g^p ) f'Ap + 4AV(2 cos k/8‘^^ + oos^kM)] ) (4.38) 

SiiTiilar expiBSsions can be obtaired for incidence along (lOO) and 
(ill) directions. Table I gives the number of partial waves contri- 
buting to the total scattering cross section for various directions 
of incidence in simple cubic, body centred cubic and face centred 
cubic lattices. 

Table I 


Crystallographic 
«xes of incidence 

Partial vJavo.s Contributing to C>-r 

Simple cubic 

B.C.C. 

F.C.C. 

(100) 

S, ? and D 

S and ? 

S. ? and 

(MO) 

S, ? and D 

S, ? and D 

S- 

(111) 

S and P 

S,?,D,F 

S, ? and D 


In the limit of long waves the dependence of the total 
scattering cross section on the incident direction ceases to exist, as 
expected, and the expressions for the two lattices are ( in units of 
a2), 


0 ^ TB 


4 , a? ^4 / AM .,2 ^ 


TT I F 


SB 


6l) or M t — y 


op 


16 




^F ““ 


"^("sFr 

1 


37r( FpgP ^-^OB V 


(4.39) 
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f 


(4.40) 
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Th'-- scatteriTi^ cross . section in bcth cases noveals Hsyl-eigli— "type 
scattering when there is no resonance. F— and F-^^ are to be 

or L r 

Gvalua'ter^ in *fclie low .'freqiiency limi't • The cliange of mass affects 
only the S wave part, anrl the cha.nge in force constants affects all 
the Partial waves. 


As discussed in the last chapter, this Ha7/le igh-t ype of 

scattering does not hold at resonances. At a certain frequency 

the total cross section will be proportional to CO absence 

r 

of resonance but at resonance it will be much larger and xd.ll bo 

p 

proportionel to Cc> “ . This is because when there is a S-type of 

r 

.resonance, the Be g p will vanish and Im Fsg p is proportional 

so that ^ m is proDoxtional toT --t- 1 

iB,F oj 

Thus the total cross section becomes larger by a factor f 

^ U) 0B,F/ 

The same is true for the ?— wave resonances where He F.,^ „ 


to 


( 1 

^^0 0B,F J 


-2 


0 


and Im Fpg ^ is prooorticnal to 


cO 




^ 0B,F/ 
by eqi 

easily that |1 is proportional to ^ 


The width of 


these resonances will be gi'-en by equation (2.43) . It is seen 


\ 3 

^ 0B,F/ 

The resonances in the scattering cross- section can be 


investigated in the same way as in chapter Til. The figures 7 & 8 
show the results for body centred cubic lattice. Here denotes 
the value of Eg ( defined in Appendix I ) at resonance . 

From the figures one can see that the ? , B ard F- like 
resonances occur only when the force constants are decreased 
sufficiently. These of course do not involve the change in 
mass. The S-wave resonances depend on AM and iiV- For 
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slight incrsase in ir:£ss snc? sufficisnl ?58crs£.s^. in ths force 
constant eno obtains resonances* For vniy high value of the 
increase in mass, resonances eopear near the banc^ edge, no matter 
whether the force constant is incrsasec^ of decreaser'^. whereas for 
slight changes in mass the change in ferree constant plays a 
dominant role* The localized moclos can be investigated using these 
resonance denominators outside the band. * 

^ i^ 4_Zhpnon Scattering and the mean free path 

The mean free path of phonons can be computed using 
eqTiation (2*48) . This involves (5r^ ( 0 ) and. the scattering 

angle G » (5^^^ (O ) is readily computed from the scattering 

amplitudes. In the long wavelength lirit only the S and. ? terms are 
significant upto the order • The expressions for the 

scattering amplitudes for the two lattices in this limit are as 


follows . 

(a) Body Centred Cubic lattice: 
f 


( 4 ^ (^) 0 

. Ai 


V 


cos 0 


AX') _ 

WttJ 


V J?pg / 


(4.41) 


which, when substitute'^ in equation (2.48) , gives us, 



(b) Pace Centred Cubic lattice: 
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anf^ Procee'^ing in the saine way we obtain. 



It may be noted that in the expression for the mean free path 
the two partial wayes lead to an interference term like the one 
gxyen by equation (3.34) for the simple cubic lattice. 

In conclusion, it nay be stated that even in the somewhat 
unrealistic inod,els of cubic lattices stud.ied here, some important 
features of the scattering process of phonons from a substitutional 
impurity becorre evident. The total scattering cross section depends 
on the direction of incidence of the phonons, as also the number of 
Partial waves that contribute to the cross section. But in the limit 
of long waves, no such dependence is there due to the isotropy of ' 
the constant frequency surfaces. The dispersion laws of the models 
considerdd here lead to Debye model in this limit. The mean free 
path for long wave phonons contains terms that arise due to' the 
interference of S and ? partial waves. A comparison of the results 
obtained in this chapter with those obtained in the last chapter 
shows a great deal of similarliy among the three types of cubic 
lattice s . 



CHAPTER Y 


A MOREL ’*nTH PHONOF POLARTZATIOF 
'? 5 .1 Introduction 

The proces? of the scattering of ohonons from a substitutional 

impurity considered in this chapter includes the polarization of 

phonons in the lattice. This problem has been discussed briefly by 

45 

Thoma and Ludwig for the isotope defect. Maradudin has treated 
the isotope defect case in soite detail a*nd derived the scattering cross 
section in the limit of long waves by using Debye approximation. Here 
we treat this isotope defect case in detail without using Debye 
apprcKimation in the limit of long waves and then we discuss the 
problem for a general substitutional impurity with grouo theoretical 
partial wave analysis which has been done for the scalar models in the 
last two chapters. 

The problem is mathematically involved in the sense that 
significant part of the calculation requires cumbersome numerical woik. 
Physically the process is conplex because here the concept of polarization 
of phonons is not as simple as in the familiar case of photons. In 
general each state of polarization of phonons corresponds to a 
dispersion law peculiar to that state, unlike the case of photons 
where all the polarization states correspond to the same dispersion. 

All these dispersions become identical if for example the model is 
such that one Cartesian component of the displacement of any atom 
is not influenced by the other Cartesian components of displacements 
of the same atom or its nei^bours. In chapter III we considered. 

such a model. 
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The scattering problem can be stated as follows. Let phonons 
with a definite polarization and hence corresponding to a dafinite 
frequency branch be incident on a defect site and get scattered. In 
general they will get scattered into other states of polarization 
also. As expected the scattering process depends strongly upon the 
incident direction. It is required to find out the scattering cross 
sections for the process. In the model that will be considered here, 
the constant frequency surface at low frequencies will not be spherical 
as is the case with the usual assumption of Debye approximation in 
this limit . 

S 5 .2 The model 

J 

A monoatomic body centred cubic lattice is considered in which 
each atom has mass M and is coupled to its nearest neighbours 
through central and non-central harmonic force constants ^ and V 2 
respectively. Such a model has been discussed by Fine and subsequently 
by Launay^^ . The main assumptions and approximations involved are 
the same as discussed in chapter II and the dartesian Co-ordinate 
system is identical to that described in the last chapter for the 
body centred cubic lattice . A general substitutional inpurity is put 
at the origin whose mass is M + ^ M and which is coupled to its 
nearest neighbours through force constants ■Vi*''AV]^andV 2 
4-AV2 ' "total number of atoms is taken to be N. and ^ 

are supposed to be unequal because if they are equal on© gets 
results identical to those obtained in the last chapter. This is due 
to the fact that when ■Vl ' ^2 , then for example the Cartesian 

X - component of the displacement of an atomr is not coupled to 
the Y or Z components of the displacements of the same atom or 
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its nearest neighbours. Thus to make a transition to the scalar model 

from this, one has to set = V 5 . 

The time independent equations of motion can he written in 
47 

the matrix form ^ ( ve shall write ‘is- U {E ) / in this chapter) 

DU - M ^ IJ = 0 ( 5 .1 ) 

A 

where D is 3N x 3M matrix and U is 3N - dimensional vector. 
shall follow the following partitioning. U is partitioned into three 
U - dimensional vectors = 1,2,3 such that has the 

X-displa cements of all the N atoms as its elements. Also D is 
partitioned into the corresponding N x N blocks D^j^ each of which 
is a cyclic matrix. Using £ , m, n defined through equation (4.1a), 


1 .e . , 


( £ + -g — ) ® “g" ® ® 2 ^ + n) a e 


and defining the two constants 

Vi + s V2 




the elements of 

D^|j can be written as follows. 



ii 

r ^ 

^ ^ m\m 

Sn',n 

^ S ^ IJ \ 



, 1 0 m' ,m + ij 


*'■ S m' ,m - 1 

1 

* £-1 SmSm + 1 

S n* ,n - 

for oC = 1 

f2 and 3. 



.m+l 


,ni+l 

( 5.2 ) 
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D^2 ( £ Ein, £ 'ni'n’ ) = ( £ icn, £'m'n' ) 


21 


^^2! ^£',£ S n' ,n n'm + 1 6 m+1 (S' n]n-l 

^ m'm -1 n'n + ij' ^n'n^ f+ 1 m’ m - 1 

fc ” 

^ l~l 1^ m^m + 1 1* / ^£’ £+ 1 S'm'^ m-l S n'^ n + 1 


. £-1 O m’m + 1 <6 n'n - 1 

> ^ i j 



(5.3 ) 


13^2 ( £nin,£'in'n' ) = D ( /mn, £'m'n' ) 


sj ^ njn ^ m'^m+l £’ m'm+1 Sn'^n-l 

^ ^ L' L‘~^ Smjo+l 5 n'^n-1 ^ *" Sm’mr-l 

^ ^ U £“^ ^ m^m+l I 


and 


(5.4 ) 


^23 ^ ^ £'ni'n' ) 


^32 ^ X'm'n’ ) 


/^2j ^£'J. ■5” n^n-| ^ m’^m+1 '*’ S m'^m-^ " m'm+1 S'n^n-l 

6^ m',m-l S iii',n+l^ (l\ £+1 S m^m-1 t> n’,n+l 

"■ S i\ £-1 ^ m'm+l <fn '^n-1 I ^ 5n^n|S£* £+1 5'n_;m-l 


^ ^ li\ JL~^ S^mjm +1 


(5.5 ) 


The soluticsn of equation (5.1 ) is denoted ty 2o * 
%bstituting i 

(R) = exp (i k - . R ) 
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in equation (5.1) leads to the equations 


whe re 


(k ) 




M <j0 ^ 




( 5.6 ) 


0 = 2_i /2 (E ^ S' ) exp l± k.CSS - S )? 

oLf^ o(-p U J 


defines the elements of the dynamical matrix (T . This equation 
(5.6 ) becom.es , 




M o..r 


V 


(5.7 ) 


unich is an eigenvalue equa.tion. This gives three eigenvalues 


MCA3''(k),3 = lj2,3 which are the three frequency branches and 

J 

the corresponding eigenvectors v denoted by e^^'^(k ) with 
components (k ), 2 tcC - 1^2 and 3 ; are the polarization vectors 

o6 

for phonon propagating with wave vector k . For the model under 
consideration, the Cj” matrix has the explicit form, 


(T = 


^ ( 1- ) 8 

8Pi(l-C^CpC2) ^p2^1^2'"3 


^ /3 2 ®i^2^3 


P/3p. 






(5.8 ) 


with = Fin ~ — and , = cos -g — 

O 

The eigenvalues M CO ^ (k ) are found from the roots of the cubic 
e quati on , 

At- + C^s|s|A. 2 ( 8 ^)=(s 2 s=s|c^Cj,C 3 )= 0 ( 5 . 9 ) 


with the roots 7(.? •? = 1>2,3 related to the eigenvalues 

J 

Aj = 8 Pi ( 1 - C.CpCj ) - M (jO ^ (k ) 
irhen the eigenvectors -(k) are 


(5.10) 



3 0 ')(k)= 


'kfiz (®i 


¥3 " hhh - ^^ 3 ) - ^ 
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8^2 C,S„S,1 /d, 


L 


a (=i=2=^3 * \¥3 - - »^a 

2^3 ^ ^lV3 - ^^3) -Aj}[A. - 8 ^ 3 ^C^}/d 


with j 1,2,3 and D. aiB proper nonnali 2 ation factors. 

j 

ihe polarization vectors satisfy the orthonormality and closure 
islations , 


2 

(5.11) 




end 


z .-jS’ ■ s.f 


( 5.12) 
( 5.13) 


( i ) i) 

and products like e^ (k) e^ (k ) transform as k^ k^ under the 

symmetiy operations of the point group of the ciystal under viiich 
2 

CO (k ) remain invariant. 


’'’'ow let the substitutional impurity he put at the origin of the 
co-ordinate system. The equation of motion will now read as, 

DU -MLO^U = PU (5.14) 

whose general solution will be 


A A 

U = U. - GPU 


(5.15) 


where the Green function G ha:s‘ elements ■ given by^'^ 

G; (og^'b b' 1 =i_y ri! 

dA ^ -~W? (’kr+ i7~ 

^ ' ( 5.16) 

2 

In particular^ ^ cliagonal in 06 , p for 

the present mod.el due to the transformation properties of 
®^*^(k) e^^^(k) under the qymmetiy operations of 0 ^^, the point 
group of the ciystal under coftside ration. 
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Thus 


47 





^ “(jJ . (k) + i 6 

tl 


(5.17) 


vhich shows the equality of the three diagonal elements as well. 

5.5 The Isotope hefect : 

From equation (5.14) cue gets 

u . (R) = u (R) - / G , P>3.-,/(R',R’’)?)a^(R’’) 

^ (5.18) 

Now let the substitutional impurity be an isotope of the noimal 
atom. Then = AYi = 0, and the pertuAation matrix P has 

non - vanishing elements corresponding to the defect site ^ich is 

taken to be the origin. Noting that G and ? have been partitioned in 

A 

the same way as D one kas ^ 

P/S/3' (a'.B" ) =AMtj2S/3/2,' 6 'e,,o Sr,o 

Substituting this in equation (5.18) one gets 

U^(R) =noo( (S) - G^.(ca? R,0) (0) (5.20) 

A 

Putting R = 0 in this ^ one gets ^ 


''cC (2) = 


■:v- 


1 + Am to 2 G (U) ^0,0 ) 

2 

where the fact that (CO ^0 ,0 ) is diagonal in oC and /S has 

been used. Substituting equation (5.21) back in equation (5.20) yields 
the result . 


( 5 . 21 ) 


- 


” AM oo 
j i + AMOO^ G( C0^0,0) 


o.^p(u) ,S.fi) (£) 

A (6.22) 
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As discussed in Appendix IT , G^.-^ ( U) , R , 0 ) can be expanded 
for large R as 

^0) = -2j ^ ^ exp (i k^. R ) (5.23) 


where each term coriesponds to a frequency branch denoted by the 
index j . Substituting this eouatibn in equation (5.22) one gets 


U , (R) = U, 


^ exp (ik^. R ) 


(5 .24) 


whe re f ^ = 


Amo:^ 

1 + AMt0 2G(o3f 0,0) 


-- 2 -^^(2) 


(5.25) 


‘iere the index i refers to the polarization of the incident phonon, 
'-^hus we have the solution for the scattering problem in (5.24) for 
the isotope case. ,.0ae can fix the incident wave by choosing _TJ^ 
suitably. It is clear that for a given incident phonon with a given 
frequency and polarization, the scattered wave consists of in general 
three polarizations. Thus an incident vibration propagating in a 
certain direction with a given frequency and polarization is scattered 
from the substitutional impurity and at a large distance R the 
resultant vibrations can be analized into three distinct waves 
trave-lling in general in three different directions with three distinct 
dispersions corresponding to the given frequency of the incident 


vibration i,e with three different group velocities. 

p 00 


OO 




Incident 



- to 


, Ui V tO , 


3 k3 

“o 


Impurity 
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We shall call the three final Polarizations as the three 
polarization channels. Equation (5.24) permits one to separate the 
Contributions to different polarization channels. Thus , if one is 
interested in the scattering to the first channel ( with j = 1 ) only 
then this process is described by 




' Wi R 


(5.26) 


Similarly for the channel j = 2 , one gets 


U^(R) - ^ exp (ij^.H ) (k,k^ ) (5.27) 

NX R 

and so on. Then these scattering amplitudes can be used to find 
contribution to the total scattering cross section and mean free path from 
phonon scattering into different channels. 


However one can employ the optical theorem to find the total 
cross section <ri directly. This involves the imaginary part of the 
forward scattering ' a-rplitude ■ i which for the present case becomes. 


d (k,k) = 
3 


AM tO 


1 + A-M6d'=^G( U) 70 , 0 , 


) e*'J^(k ) 

P P 


(5.28) 


where k is the incident wave vector and j is the incident 
polarization index. 

In Particular , let the phonons incident along (lOO) 
direction be considered. Then k = (k, 0,0 ) . Further let them 

be polarized, say lengitudinally with frequency given by 64) ^(k) • 
Then we have 




(5.29) 
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and 




/ A M o) ^ 

\ 1 + AMCd%(LO ^,0,0) 


) 


( 5.30) 


In the limit of long waves one gets for this incidence, 

y?“' • >« • > >] 


= ilL 

k ~ 


AM u) ^ g^(k) 


U' 


(l+ AMtO^G( LO ,0,0) I J 


r 


2TT^ 

But g (k) = — 

1 /3,v* 


167T 

M V*U) ^ 


- AMoO^ Im G(U> ?0,0 ) 


(5.31) 

(5.32) 


where 00 ^ = 8 ^i/M 


Also Im G(t>0 , S ,0) = 


( 5.33) 



TT to 

6 M Y* (O ^ 


for small to 


(5.34) 


with F 



X cos ( -^ . ~Yi_ ) ^ 


(5.36) 
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^ (27)2 sin Q cos^© cos^w sir? (p 

where cos = L. I 

(cos^© + sin^e cos^q^sir? Cj? )^/2 

We have also to use the feet that for , 

k = (k,0,0) , M(j0^=ygk2 orCO^ = 

/I 8 


2.2 
s’ 


so that 


k = - 


(5.37) 


Thus finally ve get 




4 ^ A M 2 / u3 ^4 °l 

TT M ^ cAo^|1+ AMU)2g(ID2,o ,0)j ^ 

in units of a^, 


where 


4vaTr^ “] 

b = 

r rr 7 

L 6V^2 J 

1 

_ 24 X 16 TT j 


(5.39) 


The San® result is obtained for all initial polarizations with 
(loo) , (OlO) and (OOl) incidences. This result goes over to that for 
the scalar model given by equation (4.39) of chapter IV ( one has to 
set AV= 0 in this equaticxi to get the result for isotooe defect 
case) if we set 2 "0 = y 2 so that ^ ^ “ 127fx 16 \ffi 

from equations (5.35), (5.36) and c^^ is unity. This transition from 
the present model to the scalar model by adjusting the force constants 
suitably is interesting to note. 

The expression for above shows "Rayleigh type of 

scattering with possibilities of resonance at idaich LO-2 

The remarkable point that emerges from this is that although the 
dispersion used here does not go over to Dehy» model dispersion in 
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tht, limit of long waves like the scalar model, we still get the same 
result except for a numerical factor. Thus our result differs from that 
obtained ^5 using Debye model by a numerical factor. However, the 
difference lies in the fact that this numerical factor depends upon the 
direction of incidence o.'aA the polarization. For example , for (llO) 
direction of incidence , c^ in equation (5.39) must be replaced by 

0 1 in a typical initial 

choice. This refrects the deoa:.i;ure of the present model from the 
Debye model in the limit of longwave s 

^ -5»4 General Defect : 


Now we consider a substitutional impurity with changes in mass 
and force constants. Let M + A M be the mass of the impurity atom 
which is coupled to its nearest neighbours through force constants 


be partitionec^ in 

the same way 

A 

as G 


j — A 

A 

A - 

A 

P, 

P 

P 

p = 

11 

12 

13 


A 

A 

A 


P 

P 

P 

1 

21 

22 

23 

1 

A 

n 

A 

A 

1 

_^31 

^32 

^33 _ 

and each of the 

N X N 

A 

matrices P. 

partitions p. . in the 

r- 

space j 

sites. Then one has 

A A A 

A 

A 

^11 ^22 

and 

66 

p. . 

“ ^ - 

^ • 

with p^^ (O ,0 ) 

= iiiM - 

■ 8Ap 

= Pll(Sn,S) = 

APi, 

Pll (Bn »§n> 


A 




ij 


"11 


(Q »Q. ) “ Pji (Bn »5n ) ” 0 , i / j 


( 5 . 40 ) 
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ana (0 , Sn ) = (Sn ,0 ) = + A j3 2 . i 3 

where the sign depends on i and j and the nearest neighbcsur 
whose lattice site is given by Rj^ . 

He re A p ^ = and A p ^ = 


(5.41) 


We follow the partitioning already described and define 

A A 

S=PU,Sq=PUq to get the scattering equation as 


A A 

S = So -P G S 


^ck ~ 




whence s - 2^ p^^ 


( 5 . 42 ) 

( 5 . 43 ) 


where s^ is the only non-vanishing part of and lies 

in the space. 

If we define a 27 - dimensional vector , 


s = 


and 


A 

g = 


^3 


A 

P 


A 

A 

A 


^12 

^13 

A 

A 

A 

^21 

P22 

P23 

A 

A 

A 

^31 

P32 

P33 


A 

gii 

A 

%2 

^13 

A 

A 

A 

^21 

®22. 

^23 

A 

A 

A 


%2 

®33 


then w© get 


A A 

S = Jo “ P S ® 


:( 5.44 ) 


whose solution is 


it \-l 

s = (I + Pg ) Jo 


( 5.45 ) 
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^ 4- A A 

= V M V 


(5.46 ) 


«here the 7 - mtrtx found group theoretioaUy block diagonolioee 
(I + Pg to M. 

chh Qs in. ciisptor XI. writ© M £><? m w 

J ap ei. ij., wnxe JM as h - Z_j corresponding 

to diffeient synmetiy blocks and get 


so that 


A , A A 

£a7 = V V s ^ 





From these we obtain expressions for 
known and then 


-V 


(5.47 ) 


whence S ^ are 



(5.48 ) 


gives 


^oC (S) ) S^(Rn) (5.49) . 

For laige R one can expand ^ ,R , ) asymptoticallj 

to get. 


V(R) - U 





j 



(5,50) 


where k and i specify the wave vector and polarization of the 
incident phonons and F (k^ kj ) is sons trigonometric function of 

k and ^ . Then one can use the optical theorem for the total 

scattering cross - section. Iw addition to the asymptotic expansion 
of Green function , one needs the elements g ^^ ( ^>Rn#Rn ) 

and the V - matrix for explicit evaluation of the scattering 
amplitudes F (k^, ^ ) and their partial wave expansion' . life 



representation 


proceed to discuss these in the following. 

To find the V _ aatrix we mst find the 
r of 0^ in the space of the displacements of the impaiity a„d its 
nearest nei«hhcnrs. For the model under ccnsideiation , this space 
is twentyseren dimensional. Ho*yer it is the di„ct product of a throe 
dimensional space consisting of the Cartesian co^onents of a displacement 
and a nine dimensional space consisting of the position co-onlinates 
of the affected lattice sites. Hence T is a direct product of T ^ , 

the representation of 0, in the ( x . y , a ) apace am ’ 

Is the same as that for the oorrespendi,^ scalar model uhich for the 

present case is =21 + 1 ? 4 .ir mu 

2 Ig lu 2g ■^2u • Then we get the 

result also obtained elsewhere in a different way, 

(2 At_ + F-, + P„ + b \ 

Ig lu ^2g ^'2vi ^ 


= F 


lu 


'‘ig " ”lu * ®g * 2 "'ag * *2u * ''ar + % + En (5.51) 


where we have used the results , 

^lu®I'lu ‘ '‘ig * \ " ^Ig + ^2g 

J'lu® fsg = ■‘su + Eu + 

F. 


( 5 . 52 ) 


and 


2g 


The above equation inmiediately show's that in the partial 
wave expansion, we will have S, , F^, G . and H type 

of waves where the nom-ehclature is that depicted in Appendix III. 
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The sum of and X >2 waves in the R wave and that of and 

i ?2 is the F wave . We have distinguished between and II 2 because 
they arise due to different irreducible representations of . 

The Y - matrix can be evaluated by using the projection 
technique des>*.ribed in Appendix III with the aid of equation (5*51) 
for P , An alternative method v;ill be to use the C'l.ebsch 
Gordei Co - efficients for the direct products given by equations 


(5.52 ) and use the known basis functions for P' 

A 

columns of matrix 7g given in Appendix III alcng 


the basis functions for 



simply 



which form the 

2 

lidth the fact that 



This method is elegant and useful. The resulting matrix is depicted 
as Vp in Appendix III , the columns of which are eigenvectors of the 
symmetry types in the following order , 


The rows are labelled in the same way as those for the 

A X 

matrix p . In fact Vp can be partitioned in the same way as p to 

give , 



and a conparison with Vp depicted in Appendix III gives the 

A 

structure of these 9x9 submatrices V^^ . 


For example , 
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0 

0 

0 

d 

0 

0 

d 

-d 

0 

0 

-d 

d 

0 

0 

-d 

-d 

0 

0 

d 

d 

0 

0 

d 

-d 

0 

0 

-d 

d 

0 

0 

-d 

-d 

0 

0 

d 


with 


1 

7t 


1 


(5.54 ) 


(5.55) 


oimilnr ly one can write down other submatrices. 

It has already been pointed out in chapter II that the 

stmeturo of the block diagonal form of p , g , (l + Pg ) and 

, .-1 

{ I ■+ pg ) can be known from the structure of the reducible represen- 
tation P . Thus from the expression for P given by equation 
(5.51 ) and the order of columns of the matrix 7p described, we 
expect the following block structure given in order. 


Three equal 3x3 blocks for P- wave. 
One single element for S - wave . 

Two (iqual single elements for - wave . 
Thrno equal 2x2 blocks for I >2 - wave. 
One single elisment for " wave. 

Three equal single elements for F - wave. 
Throe equal single elements for G - wave . 


Two equal singi^® elements for H - wave, 


The structure of Sn > ISn ) written 

down using equation (5.15) and noting that (k ) ) 


transform like k^ 
Also wo note that 



under inversion and greup operations. 

* f 

(k) (k) is real and If we set 

- P “ 
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Sn = Rn-M , ^ \^\) 

This leads to the conoluoion, (we suppress cO ^ xn t ie foll'~^wing) : 

g^|^(en > (Rn , ) (5.56 ) 

^^(ptSn . Ei ) = ^ (ai , Sn ) (5.57 ) 

For the model under consideration, 1 Rm | have five values 

which appear as subscripts 0 to 4 in the following which in terms 
of the components read ( C , 0 ,0) , ( + i , t 7 , 1 i ) , ( ± 1, 0,0) 
and its cyclic permutations, (i. 1 , 1 , 1 ^ 0 ) and its cyclic permutations 
and ( + 1 , + 1 , + 1 ). These are the sans as those for the scalar 
model. Inserting these values of into the expressions for 

2 L- 

('*^;|Rnil ) we get the following different Green function 

elements, 

,RA ) . o( = 1, 2, 3 we have 

to ^ _ oO? (k ) - ie 

3 “ 



3 N M 




LD ^ (k ) 

3 


( 5.58 ) 


g- 


D 1 


N M 



e^^(k) cos ^k^ cos cos -^kg 


a> 2 - (jO? (k ) - i e 




k 3 


cos cos i-kg cos - |kg 
(jL)^_U)^(k) - i£ 

J " 


3 N M 


(5.59) 
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^ ST* X*" ^ ^2 ^2 

“" ~[3-r7-j)T(7~:7g 

£. 3 D 


1 ^ cos k^ cos kg cos k^ 

3 N M — I 

k 3 3 “ 

„ = J: ^T V V V 

D21 N M £^ L^ 2"" J’a ( k)'!'! e" 

k j 3 


0^ = 


1,2,3 


g. 


•a-/" ■; ') ( T ■) 

1 -ST" ^ (k) (k) ^ ja 


2.2: 


D22 N M Z^ Z-* (Ji)2 _ (k) - i 6 

2 j j 


g. 


oC f ^ ~ ^» 2,3 . ^ ^ oC 

*(i) ( 3' ) 

e (k) cos k cos k^ 


D31 



for (/ j ^ = 1>2,3, o^-. / p. 


and 


’D32 




a*(3)(k) e:j)(lc) cos k^ cos k 




(4? 2_ 2 ^ g 

cL f ^ f ^ ~ ^» 2,3 

such that 

and d ^ • 


( 5.60) 


(5.61 ) 


(5,62 ) 


(5.63) 


(5.64) 



* • 


• « 


ns 


The 


are the independent elements cccuring in ) , 

^22 (Sn ^ Rfi ) and (Rn » ) . 

j R^ ) , 0 < ^ ^ ^ ^ = 1,2,3, we get 


%1 - ^ 

j 


1 X" (k) e t 7^2 ii-kj) 


CO ^ ^ 


g- 


N4 


N M 


II-' 


- a)" (k) - i £ 


®{i ^(k) e"-^- ^1 i ^2 ± ^3 ) 


OO 2 


- U) 2(k) _ i £ 

D 


(5-65 ) 
(5.66 ) 


®N3 - 


TrZZ 


S ) e<i> 


(y ±^p ) 




U) 


u) : (k) - i € 
J 


(5.67 ) 


The subscripts D and N infer to diagonal and non-diagonal 
Partitions of g- The first index refers to the | | values and. the 

second index refers to the different types arising from the same 
ISml value . 

With, these one can block diagonalize g and p and hence 

. aa . 

N Pg ) vdiose xnverse can be found and then s f vectors can 
bo calculated . 

The resonance denominators are given by. 


F. 


r 

^ A ^2 ** ®D4 ^ ~ ^^D21 ” ^D3l ^ ^D22 ~ ^D32^ 



(5.68 ) 
(5.69 ) 
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where 




8cA)3^ 16 A^. 

' ^A'l - A Pi 0 


i6A^^a 


-AS 


16 d g 


°°^D1 [gpo ^ ^154 + §221 ^ ^D3l J^(gN3 gN4 ^ 

+ 2(gD22 ^ &D32^'3 

i 

ISdgjjl Jl( gjj3 + 

“ SN3 %4 "I 


where c and d 


are defined in equation (5 .55 ) 


1 - A^i j (§20 - § 24 ) - (§221 - gD3i) + 2 (§ 22 . 

^ ^ %4 M 

y— —J 


D22 “ ^D32 ^ 


(5.70) 


^^2 1 ■ ^D4^ ■*■ ^^D21 “ “ 'SjJS ^ %4^ 


kJ ~1 (g 22 - §24) - (§221 - §031) - 2(§222 “ ^ 


2( A 1^1) ( gjc “ gN4 ) 


( 5.71 ) 


Fpi - 1 - (A^i ) |^(§ 2 o + § 24 ) + (§221 + gD3l) - 2 (gD 22 ^ ^^ 32 ) 

" - %4 ) J (5.72 ) 

^F2 ~ ^ ^ ®D4^ “ (^D21 ^D31^ (®N3 “ ^N4^ (5.73) 



Fg = 1 - (A^j) 


(Sdo ■ *D4> * - Smi' * *^^113 * ®N4^ 


(5.74) 


and 


J'h = 1 


- (Ap,^)^gB0 ♦ S^i) * (BdEI * «wi^ ■ ^ 


6d32^ 


(%3 ” %4 


) 


(5-75 ) 


An alternative way of finding the vectors xs to use 

the _ Gonion Co-efficients and to make use of the results for 

the scalar model. This is easy to do with the single factor blocks 
for bigger blocks the problem of indexing the Green function makes the 

t*) po s 3 oinoXi C3 'fcGci- • 

once the vectors ere calculetea , one cen find 

F (ki , !d ) aopearr^g in emotion (5.50) . The expression for 

F (P , kb oontsins contributions from S , p, , D^, F^, Fg , 

■, H„,^ver in the limit of long uaves, 

0 end H typo of waves in general. Ho-. 

the F ,F ena H type of waves do net cont^bhte . This can be 

^ ■ ■ the result for the scalar moJel where in the 

anticipated by examining 

nlir the waves corresponding to 2 k-j^g 

limit of long waves only the wa e .^^tions 

. w nee in the present case one expects to find ccntrrbutre 
contri-butn . Hence in me pre 

from waves oorrrrsponding to F^^ ® (2 %g * V 


= 2r,„tAigtFgtF^tF2g 


(5.76 ) 

i '.J 


■ T those waves that contribute to F(k » k^ 
and this leads to preexsely those 

tn the limit of long waves. Thus in this limit 

i 1 a F^^bkbk^ ) ^ ^ ^ 

F (k"-,ki ) = ^ ^k ; 

fn.\ i 4 .y 


(5.77 ) 
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where k and i specify the wave vector and polarization of the 
incident phonons and and j specify those for the scattered 

phonons when k^ lie along synmietry directions • Using the long 
wavelength limit expression for incidence along (lOC)' direction, 


g ^ (d) = 


_ 

"id : 

TT X 8 ^ P 

T? M u)^ 

0 


with (jQ ^ 
0 

M 

1 

tho exprepsions for 


) are found to 

-Jf; 

F^^d . = 

ki ) 

/ 1 ) 

/Ah)^ 

y 

\ 3 TT Fg ./ 

\ /sj 

L 

L a 


( 5.78 ) 


/ - 0 t 
^ I' 


(5.79) 


''‘’HkSkj ) - (hi + Sb^)[yelf'(k) )] 


Tr mlO^ Fp 


I 


(5.80) 


with 


h^ = A M jAgg A33 - A23 ^32] ^ 

-^|^^12A23 “'^22^1^ 

A 23 A ,31 " ^21^^ ^ Ap 
X 1^21 A 13 - A 23 A 


(5.81) 




Am lO 


(5.82) 


' ’ A V A . Jv A 

where A= I + Aji , the 3x3 matrices ZXp and o 

Jt O 

are defined after equation (5.69 ). 


g 
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whe re 


^ 4 ~ ^ ^D21 “ ~ ^D32 ^ 

(5.85) 
X 



In all the equations above one mst evaluate the quantities 




2TTFJ 


A^l 


kgi^fk) - kjigtk) 
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■'S ' ^ ^Dl > ^-02 ’ weve length limit. 

Let the fincideny phonons be incident along the (lOO) direction. 

Then k. = k - o i 

^ '2 j kj - 0 and the total scattering 

is given by the optical theoiem. 


cross section 


CT- = JTT 


Im 


F(k\ k^rj 


(5.87) 


(t) j^ngi tudinal Polarization Here (k) = S &i > 

p — r 

I'ow F (k , k ) with k and i specified as above leads to the 
following partial wave contributions^ 


, kO = , 

3TTFs 

2 k^ 


,(P) 


Tf Mu>% 

O ^ 


(hi- 8hp 


3Tr F, 




= 0 


F 


D1 

(G) 


(5.88) 


(ii) Transverse polarisation Here e^^(k) = 5^2 or SliZ 

f(s) = 0 = f (dO 


.(P) = 


7T Fp 


(hi +8h2 ) 


" ^ ‘ - ri-;- ' ' 


(fi) 

F^ ^ = ( - 


2^^ 
1 


D2 


(5.89) 


) ( ) k2 

2 71 Fq f9>i 


Now we ha VO for the Green fanctions in the limit of long waves 


JCj bj 


®DO ~ n ’ ^la. “ 48 ’ ®D4 Kj"+- 


"l2 


(5.90) 
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with b 


?. 


z. 


r 


-j 


sin Q Q d (p 

r ’f) T ^^2 


( 5.91 ) 


horo (e ,(p ) is given by equation (5.35), 

Tm + x^b^ ) , im ) 

%?>2 " ^1 ^2^4) » gj)32 = ( + Xgbg ) ( 5 .92) 

Im =-■■ Xg , Im = 4x2 ^7 > %3 = ^2 ^8 

where x^ = - ^ Tf 

^ 6Mv^a}^ ’ ~ 

0 ^ o 


and , i = 3 to 8^ are cumbersome integrals over angles like 

b^ and bg described inequations (5.35) and (5.91). These jaeld 
numerical values after integration with given ratios of ( ). 

/3 1 


Using all these results and the optical theorem along with 
the fact that for (lOO ) direction of incidence, = ~Z_ ^ one 


!A)‘ 


tho rosiilts, 




92 


* c 


4 - 




&.rj 


h9 2_^ 

fzl 


2 



/Afi2 \ 

I 1 ^ 2 / 


1 


(5.94-) 


p{ere T and L refer to "ftie transverse and longitradinal polarizations 
reapeetivcly and the mmerieal cfuantities . a^,-'i-^= 1, . . . , 10 are 

appropriate eoinbinations of i = 1 , > 8 and other ntimericci l 

factor?' . 

P irn lav expressions can be derived for other iheident 
dlToctn.ons and polarisations. However the partial waves contribatang 
for different infeident diiections and polarizations are- different. 
Table II below gives the result in the limit of long waves. 


TABLE II 


Direction of 
Inoidonco . 

Polariaation 

Normalized 

Polarization 

vectors 

(transpose) 

Partial wav 
es contrib- 
uting to 

CTt 

Partial wa"?^’ 
es interfe- 
ringin the 
exp re csion 
for Ar'(^^ 

. Bemarks 


100 

Longitudinal 

Transverse 

Transverse 

;i 0 0 ) 

(O 1 0 ) 

(O 0 1 ) 

S,P,Di 

P >1*2,0 
P>I>2>G 

S-P,P-D^ 
P_D2,P-G 
P-D2,P-8 • 

The two 

transverse 
branohe s 
are degenerate 

110 : 

Longitudinal 

Transverse 

Transrerse 

^(1 1 0) 

(0 0 1) 

J:— (1-1 O) 
.'2 

S,?,D2 

P,Di,G 

S— P ,P~B2 
P_D2,P-G 

?-Di,p-g 

Transverse 
branches are 
non-de ge ne rate 

111 

Longitudinal 

Transverse 

Transnrevse ■ 

-4s-(i 1 1) 
-1— ( 1-1 o) 

1 - 41 1-2) 

. 

S,P,D2 

P,Bi,D2,G 

P,D1,D2,G 

S-P,P-D2 

F-Si>P-D2> 

P-G 

P-Di,P-D2> 

P-G 

■Transverse 
branches are 
degenerate 

i 

! 

4 
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It should he noted that even in the limit of long waves > 
number of partial waves contiibuting to the total cross section depend 
upon the direction of incidence. This feature is a consequence 
significant departure of long wave behavior in the present mode 
the usual approximation. It is in fact related to the beha 

of the dispersion which in the limit of small k becomes. 


lO .(is) “ 
1 


. 1^1 )/ 3 j 

■■ "iT" 8 M 41? 


for 3 - 1^2,3 


with 


OOO^X 


27 ( k2 k2 l4 


( 5.55) 


( 5.96 ) 


2 ,2 s3 




This ,ho«a thst the f«<ps>«y sarta.es in thi. 

ia... «,veleh«th Xirit ate hot sphere, as ^ the ease it 

IS s,a»d, . nthoagh this aepartuie diowed itself.™ for 
esse of isotope defect in the .^erioal factor as seen alrea y, i ■ 

fa general suhrtitutlonal ioonrity <I“tte significartly. 
fiffoots the case of a generau suu 

free Path wiTl contain interfererce e 
For example the inverse mean free pa 

4-.nl >111+ ions depending upon the 
l«t..on dlfforont Partial vara con 

afraotion of incidence and the po ^ irtarferan- 

This is shown in the above ta • 

hetpeen waye. of differant -parities' . - 

^fisions for the total cross- section 

^'“^"^"“..soriythe Pi^yocontri^---- 
fhst tho ,, aaras. The lo» s»rgy 

tho change HI foIM • 5 ^ expected Bayleigh t»e 

T y,<T+Vi scattefing shows the ext>« 
and hence long wavelengt 



or scattering in absence of resonances. However the presence of 
resonance denominators include ■ the possibility of resonance scattering 
A comparison of the results obtained here with those of the scalar 
iiindcl discussed in the last chapter shows that there is a change in 
' parity’ of the waves in the two models. For example the, S — wave 
of the scalar model goes over to the P - wave in this case and so on. 
This result is expected by simple group theoretical arguments. 

c; 5.5 i'ummaiv i 

In this chapter the model considered may be claimed to be a 
step fon*ard in the approximation to the realistic case, compared to 
those considered in the last two chapters. As expected one gets sons 
results that are qualitatively the same as those for the scalar model. 
Among these are the modified Rayleigh type of scattering at low 
frequencies with resonance possibility and the interference between wave 
of different 'Parities ' in the expression for . One of the 

differences lies in the consequences related to the- validitv of Debye 
oproximation in the long wavelength limit and the shape of the constant 
i’t '' qvK^rcy surfaces in this limit. It turns out that in this limit 
i constant frequency surfaces are not spheres in general and hence 
V'obyo approximation does not hold. However it is remarkable that for 
+■>,0 iforiope case this leads to the same result as that obtained by 
\islng Debye approximation except for a numerical factor which depends 
upon the direction of incidence. For the ■ case of a general 
substitutional impurity it has the effect of making the long wavelength 
total cross - section dependent upon the incident direction unlike 
the scalar model vhere it becomes independent of incident direction. 
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discussed ..o» ore can go auer to the scplsr nidel from 
tho model urder consideration by adjusting tire force constants suitably. 

The other rsru feature here is of course the dependence of the 

scattering process on the polarisation of incident phonons. It is 

flTOe<^y apparent that introduction of n n • x- 

incroauction of polamzation iBads to a great deal 

. f rathematioal complexity in the problem. The acautering process 

quit,, imrolvod physically also. Compared to the scalar model, 

the corresponding partial uaves appear, as exreoted^to have change.i 

thoir 'parities' . ."arhaps one can safely contend that these additional 
features imst be preeent in a realistic base. 



CHAPTER 71 


CO!JCI,USIO?I 

• « 

:-he lx^yults obtained in the present study can be sujir^arized as 

The scattering cross section and the nunter of partial waves 
i.tr., I.utin ,3 to it aepend upon the direction of incidence of phonons, 
uhus orhi biting the expected non-central nature of the 'potential' due to 
^ sub;' bihulional impurity in a lattice (Chapter IV, Table I). 

Inn resonances in the partial wave amplitudes and the behavior of 
i-',c .'.altcijiig cross sections near resonances have been investigated in 
Iho results confirm those obtained by other x-jorkers. The 
ncnotoring of both the acoustic and optical phonons in a diatomic simple 
'■aide Jattice has boon studied and expressions for resonance denominators 
are obi-oined explicitly. 

The expressions for mean free paths of phonons due to scattering 
i. '.'om uul'Btit'utioi''al impurities have been obtained for small concentration 
xC ira,-nxritieo and phonons with long wavelengtlis . These expressions are 
iiown 'to contain interference terms between x^aves of different 'parities'. 
T'lcir I'xjliavior near resonances has been studied. Tlae similarity of 
;v)t!ul'..o for tho throe cubic lattices has been discxissed. 

TJio effect of polarization of phonons on the scattering process 
h.Mf.! Ivion studied by xxsing a sxutable model of body centred cubic lattice 
■n Chapter V. It has been pointed out that the results obtained for 
i'lis c-'iSG can be balden over to the case of the so called scalar models 
by nrj, appropriate choice of force constants. The Debye model is shorn 
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• 0 


■<-'0 be inadcqiaate for the 


exact description of the scattering process 


in 


the l.mit of long „av03. ao departure from the Dobie model in tele limit 
leads to the dependence of the scattering process on the dineotion of 


incidence (Chapter V, Table II). 


ntire study has been based upon the harmonic approximation 

and nearest neighbour interaction models. The effect of anliarmonicitj.- 

on the phonons in the impure crystal is rather complicated, and leads to 

50 

finite life-time for each phonon. This would affect the phonon 

'loan free path, and the transport coefficients which depend on it. The 

prt5scnt analysis is x^alid at very low temperatures wlien the anharaonic 

cffocto are small, and would be adequate to cover such phenomena as the 

kinks in thermal conductivity curves obtained by many workers^^*^^^^^"^®. 

Tho other essential limitation of the models investigated here 

.'S tliat the impurity is talcen as a substitutional one. There is enough 

51 

experJ-inontal evidence nox^r to indicate that quite often the impurity 

a, ton (or molecule) goes into interstitial positions. The sjnnmetiy of the 
impurity ’potential' in such cases is surely more complicated than xirhat 
has teen ;,iir!umod in the present analj'-sis. A proper theory of phonon 
scattering from impurities must take such situations into account. 



APPE^IDIJ T 


GHEM FJNCTIOMS FOR SC AT, at? MODELS 


The elements of uhe Green function matrix for a scalar model 
are of the form 


V 


G(a^ , R- , R ) 

J 


1 


\{( t-* ) 1 

co^(k) - o)^- i€ 


(AI-l) 


^^^llere clO (k) is the appropriate dispersion, the components of ^ and 

Rj are measured in the mits of the corresponding lattice spacings* 

Ihon k is diiiB ns ionless wave vector and the k integration is over the 

fii'^st Brill ouin zone whose volume is For the simple cubic lattice 

k integrations extend from - irto +T'mth V* = (2 7 ^)^, the unit of 

measurement being the lattice spacing a"'. 

For face centred cubic and body centred cubic lattices the 

symmetrj’’ of the constant frequency surfaces can be exploited by doing the 

48 

■’.nbogratiou throu{^ l /48 of the Brillouin zone , vdiich in the face 

centred cubic case is the part of the zone enclosed in the trihedral 

mgle defined by (001), (101) and (111) directiensj and in the body 

centred cubic case is the part enclosed within the tetrahedron bounded 

by the planes = 0, = Tf , and \d.th their 

range G defined by 0 ’'2 and 0 <1^34 TT , It is 

2 

«aso possible to combine these wedges to form a cube in the reciprocal 
space whose sides are of length 4'7r i.e, from —2 75^ to + 2Tf' such. 
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that integration over this cube gives four times the result of 
integration ever the Brilloidnzone for the body centred cubic lattice 
icid twice the result of integration over the Brillouin zone for the 
face centred cubic lattice . We choose this cube for integration 
because of the analytic convenience it offers for body centred cubic 
and face centred cubic lattices. With this choice and a change of 
variable one obtains 



(w ,4, sp = 


-f-T 



exp 

C0^^j,(2y - ife 


o 


d k 


(AI-2) 


whereas for the simple cubic lattice 


G. 



■+ir 



[±k.(R^-R.) J d\ 

( k) - cO^ - ±fc 


(Ar-3) 


The subscripts S, B and F refer to simple cubic, body centred cubic 

.'md faco centred cubic lattices respectively. The asymptotic ejq^ansion 

is obtained from equation (Al-1) for large E = I £^*-1- I , and this 

2 

is discussed in Appendix II. Since G (.UJ , depends upon 

I I ^ follows, we shall adopt the convention of writing 

the Cartesian components of R^-R-i “ units of the lattice constant 

'' 2 

'a' of tbcsG cubic lattices. Thus for example G (it) > stands 

2 

for R--R. = ae, + ae^ in body centred cubic lattice, Gj,( C<> 

— j -I * 

. , - n D _ oCo j. -t P + 4- e ') in face centred cubic lattice and 

stands for I^-R. = ^ 32 ^ -y 

V 
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so on. Here are the unit vectors along the Cartesian axes. 

The Green functions v^'inhoccur in the expressions for the 
scattering amplitudes are conveniently evaluated in terms of the 


parameters E = 3 - ^ i? _ 4 ^ 

s - 1 - — 


2 

oS 


CJ 

oB 


and E = 3 - 
F 


o 


oJ ^ 

oF 


so that the frequency bands are defined by the ranges - 3 ^E„^ 3 , 


respectively. These E , S„ 

o Jd 


ond Ep are dimensionless energy parameters , 

For the simple cubic lattice^ the Green function ele 2 :ients 
occuri^'ig ill tlic expressions for the scattering amplitudes can be 
derived from the general expression 


G(p,q,r,Eg) 


-Yl\ 3 

p d k cos pkj^cos qkgcos rk^ 


L(2Tr)^jJ 

-If 


wno 


■,'rc _ i/(M = - 1 /( 2 ^). The integral inside the 


28 


parjmthftsis can be written as . 


E -cos ki-cos ko- cos k_- i£ i 

q 1 0 — 


(AI- 4 ) 


3 (p,q,i’,Eg) = |^C(p,q,r,Eg) _ i S(p,q,r,Eg) J (AI- 5 ) 


■Do 

r 


dx cos[EgXj Jp(x) J^(x) J^(x) 


with 

and 

S(p,q,r,Eg) = \ dx sin(EgXj Jp(x) J^(x) Jj,(x) 


J 

0 


(AI- 6 ) 


(AI- 7 ) 


Here J (x) is tlie Bossel function of first kind of order p. The 

P 

numerical values of the functions C(p>c[,r^Eg) and S(p^€l^r,E^) for the 
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values of I'd n r.\ 

. ’S8 S 

' ' ^ *•«= 0. Xo^ „^es one ,t3 

0)3 


I(p,q,r,E ) 


nT^TTcO 


bs 3>r2 TT 0)3 

oS 


(n + 2^ 

vP q + j. ) 


where we have employed the 


following method. 


(AI-8) 


IP P^(k) and F2(k) 
integral 


are two functions of k, 


then we can write the 


I = 




(AI-9) 


using the formula 1 

X - i^ 

Now W3 can further write 



+ i77 



I 



(AI-IO) 


where the integration 
solution of F^(k) = o. 


IS over the surface of constant k where k is the 
All these are evaluated at the limit of long 


waves. Explicitly one has 


0(0, 0,0, Eg) = 

<0 


(Al-lla) 




oS 



S(1,0,0,EJ = 

Cp 

cP ^ 

(Al-llb) 


4Jira) ' 

oS 

3>r2Tru)^ 


oS 


0(1,1, 0,E) = 

> + 

2 

(Al-llc) 

S 


3 /27f CO^ 


oS 


0(2,0, 0,E„) = 

- ^ + 

2 fT o)^ 

(Al-lld) 

\ 7 7 f qJ 

/r,ra>^ 

3TF 


oS 
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For the body centred cubic l-itt-i no +^, ^ 

attice, the Green functions 

resonance denominators are 


who. re 


2 

A._ = - i/m “ l/SY . J and Y are Bessel functions of 

first ,nd second hinrl. II-3 quantities GBR^. and GBI. defined above have 

“ ^33 

boon cranutod and th.;i.r nur'criccl values are a?nailab?.o . Tl.e distra- 
bnbion function of the frequency (squared) for the bedj^ cen ned cubi.c 
l/'.Ttico, wlxLch equals GBI^ oxcenb for a nuuericcl factor^ has a iogarithnic 
s?n.f7.ilo.''ity in tlio n:.ddle> 0: the be id in this scalar .:del. 

n the low freoueno''’ 'lou" wrveloi. ith) litdt 01 . can i'clloi.r tr.o. 

DO ohod indicated in equations (AI- 9 ) and (AI- 10 ) to got 


In 


■B 


’‘'K Ib’ ^ 'iB ' \ <hB- I3’ ’ I 


^2B- 


in I “sr) . ^ Xb(=>ib - 


wlf-.h ^ 2 




>dB ' O-B 

['hero ; re used in the derivation of the long wavelength lini"*" of the 
soobtcr.Lng a nnlitudes . 

'"irnilarly, far the free centred cubic lattice one obtains, 

i-T 


2 


%( Cc) ,0,0,0) = - A. 


^ r 4 


oP 






*7 


0) E^-(cos cos cos k^eos kj- 




- 4 ^ ) 


GJc^',i,i,0) = g,^ = i(VoF 


Gp(cO^, 1,0,0) - %\ 3 


H-Tf 


"2F 


M 27 r) J 


-f cos k 005 k )- i-£ 

iC o 

(AI -18 ) 
(AI- 19 ) 

3 

d k cos 2k^ 

E -(cos k,cos k + cos lacosk 
V 12 23 


1 


-TI 


+ cos k cos k ) 

O J. 

(AI-20) 
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♦ r 


0 ( <x)' 
F 


= g = 1 
4F 


i-T 


L(27r)' 


d\ 


cos cos kg cos 2k^ 


-T 


E -(cos k cos k + cos k cos k 
^ 12 2 3 

+ cos k^cos k^) _ 
(AI-21) 


aiTii 


Cr ( 

P 


- ^'aF = A, 


+ T 
ftr 


U^T) 


cos 2k^ cos 2k, 


Ep-(coo k^cos k + cos k cos k^ 

iC ^ 


J 


whore 


F 


iAm lO = - lA^y . 

oF ' 3 • 


+ cos k cos k ) - 

(Al-22) 

The long w^elength Unit for these 


aiey 


In g 


oF 


^/ip’ ^ ** ^F*tP'®2P). ^ %F = ^ ^F'tF'^“2F> 


In r = 4. A (a 


■31.' - -p-lp- %f). in 3a^p, j, 


(AI-23) 


wi'Dre 


a 


IP 


^ATcO > a = 

op 2P 


/eiro) and k = 2 cO/ ^ 

OF ^ oF 


As ncntionod above , tables with numerical values of the integrals are 
avaiJ.fdxio. The numerical computations of these from the forms given 
above’! Iiappon to be quite tedious and tine consuning. Becently Mahanty^ 
has suggested tlie use of the os-thod of Fourier series expansion of the 
Gxxien ft-jiction. This simple technique is |)ased on expanding 
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x) 


n 


a^(R) sin(n7rx) 


with a (R) = ^ 
n Y'^ 


J 




3, 

k exp(ik.R) sin(n7r E(k)^ 


(AI-^4) 

(i0:-25) 


<c 

wiierG E(k) =14 0) ^ , 

(-U and X is dimensionless energy such that the band 

is alx^ays defined as 0 -v ^ i 'Tn-.r^ -i j. « / \ 

real part of G(R, x) is given by 

+00 


Re G(R, x) 



In G(R, X' ) dx« 


X’ - X 


(i0:-26) 


Ihe use of this method leads to considerable saving of 
computational, time. It may be noted that the quantities numerically 
computed do not involve mass or force constants so that these very integrals 
can also be used for Green functions encountered in the problem of 
scattering of spin waves and band electrons fron a substitutional inpuiaty. 
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iPPEKDIX II 


^IMPTOTIC EXPANSION OF GHEEN FUNCTION INTEGRALS 


A typical Green function integral can be written 


as 


R) = S (AII-1) 


whore 






(R) “ i £ 


(LIl-2} 


, is sone function of k -kidch does not involve R and ^<} .(k) 

^jjr - _ j v_/ 


Here /t (k) 

nr 

is a branch of the dispersion. The index 3 denoting the branch of the 
dispersion is to be omitted for scalar models. The integration is ever 
the first Brillouin zone and in the limit. For large R = | _E | 

the integral can be expanded by the method of stationary phase idiLch has 

, 21 9,7,45,13 

bwCn discussed by Koster and subsequently used by many authors 

TJ-o expfmsion given below is a versicn of this method suitable for our 


Using the formula 


X - ife 


■bo 


dt 


0 


c: puation (AII-2) can be written as 


^'1 A ^ » E) = - I (E, 6 ) ) 

py m 


(AII-3) 


wharo I (ll, a> ) = ij dt|jj d\ exp^-i <p (k,t) j (AlI-4) 

with the phase (k,t)~^blCk) - cd^ - k,B (iII-5) 



* • 


Wl«n cp is statiosaiy with respect to k and t me has - o 

■aid 0. These two equations determine at least one (kj , t ) 

which specify the stationary point. Using {aII- 5) one finds tit J 

— rk 


t one has 


"^0 such that 


r i 2 

o) . r) = ct) 

3 ~o 


V cd%) 

k ' 


k = k'^ 
““ -**.0 


S / t 


(AII-6) 


Now one can expand Cf? (k,t) about this stationaiy point (k^',t ) 

by putting k = kj + K and t' = t + T so that 

u 0 

f (-’ + t' £ ^ .K. + E (An-7) 

V. T ^ 


neglocting Mgher order terms 


a) 

7>^ 


, (k) 


(AII-8) 


k = r 


■a® ofA') 


Z) k • "T) k . 

1 3 


(AII_9) 


k = k 


Quo oan similarly expand around ^ bat only the^ (k^ ) term 

will give a significant contribution because the higher teirfe will 


ultimately load to terms of the order of 


idxLch will be negligible 


when R is large. Further if there are more than me (k^, t^) that 
satisfy equation (AlI-6) then a sunmatim must be dene over all of them. 
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Thus equ^.tion (ATI-4) can be inhtten as 

I (R, a ) = ^ (kb i fd-r fff A 

aOlki.t /<j° J JJJ 

— Cj 0 ' 


Tlois can be further simplified by writing the phase in matrix and vector 


notation as 


r Z:.M *0 ^ ' -s (3^ - 

i •' . 


where T" denotes the tr&ispose. 


^ A_i 

Further if we define \> = K + i-k — 'cC~J^ 

•WM> is * 

0 

tl:ien 

T '***' V-l 


z 4 + to K 


K = -i ^ A^' ^ + t voC 


^md == d p. We have 


used the fact that is symmetric which is 


obvio-us from 


its definition. Substituting this in equation (AH-IO), we 
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where the limits of p integration is extended to + Oo ^ Extending 
the limits off ■ integration similarly, one obtains 






But from equations (jo;i-. 6 ) and (AII- 8 ) one has 




Hence I (R, CJ) ) = ZL 

all 

-0 


2 

2 TT 



(AII-12) 


Putting this in equation (AII-3) one gets the final form 93 
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Now for the scalar models, the branch index j can be omitted and 
comparing equation (AH -2) with the corresponding equation of the scalar 
models described in equation (Al-1) one finds that X ,(k) is unity. 
I’hus for the simple cubic, body centred cubic and face centred cubic 
scalar models ono obtains 


C (U3, g ) 

3,B,F 



S (k ) 

Slk S,B,F ° 

0 


exp (ik^'R) 


R- 


(An-14) 


whore evaluating Q , and as defined by equations (AII-8) and 

/ i n 

(AII-9) we get 


g (k ) 




2 2 .2 
sin k + sin k + sin k 
^ 0^ 


2 2 .2 
L sin k cosk „cosk ,+sin k „cosk ,cosk .+sin k cosk cosk 
01 o2 03 c2 o3 ol o3 ol ^ 


(iII-15) 


= 


B 


1 

STY 


. 2 2 2 2 2 ^ 2 2 2 2 _ 
^2 ^3 "^^2 ^3 ^1 I 




J 


(AII-16) 


where C- = cos ik . and S. = sin . 
1 01 1 01 


For face centred cubic 




8Trr 




2 2 2 2 2 2 
\ (Cg+Cj) +Sg (Cj+C^) +S3(C3_+C2) 


r 2 2 2 

x-js^ (f’2’^^3^'^2 C^3+^1^^3 (V^2) 


(AII-I?) 


For these scalar models , not more than one k is obtained as solution 

— o 
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of equation (AII-6) in the limit of long waves because the constant 
frequency surfaces are splierical in this limit for these models. So 
the summation over is to be omitted. We have omitted this simmation 
also in other cases where we are not interested in finding out numbers 
or exact results but require only qualitative aspects of the results. 

For the case of the model with polarization 


(k) (k) where ej^^ (k) is the ^th 

component of the jth polarization vector. Then one gets 


-a) 




(j) 


^ (kj) (ky /R 

j all k^ ' 


(Alibis) 


wliere 



(AII-19) 


V 

The vector O and the matrix OO in equation (AlI-19) have to be 
evaluated at k^ . 


• • » 



appendix III 


EViLlTATION OF 7 - MATRICES 


Consider e lattice whose symmetiy point group consists of elements 
with the corresponding operations 0^^ . Let N be the order of 
the group and let there be h irreducible representations \ >V= 

of this group, tfe denote the dimensionality of the irreducible i^epresent— 
ation I ' ^ by n^; and the unitaiy matrices corresponding to operations 

Oj^ in this irreducible representation by ^ ^ (R) whose traces or 
characters are given by Hr) . In addition to this point group 

of symmetiy, the lattice has translational invariance but this is 
destityed when we put in a substitutional impurity and we assume that 
the ^stem still has this point group symmetiy. ,The displacements of the 
substitutional impurity and its neighbours form a space S Tdiich is 
invariant under the group of unitaiy operators • The characters 

")C_ (R) of this space S with resoect to the group can be found and 
t5ien. the representation pi of the group in the- space S can be 
written as , 




(j>) 


where 


42 


V 


- 




N. 


g 


( A III -1 ) 
(A III - 2) 


R 


This shows that the space S can be split into irreducible 
aubspaces S ^ ^ ^ which occUr times for a given y . If there 

is a non-zero vector r^in the space S then the operators, 

(A III-3) 


p ( ) = 11 z: d' ^ ’*(R) 0. 

rs N rs 


R 


pro.iect out those vectors in the space S that belong to the 
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row of the V irreducible representation of the groupf^ This is 
subject to the condition that ^ ^ f. is not identically equal to 


X. 2 


A( ) 

Since D ere x matrices, there are in all 


n . = N such projection operators. Using these we can form the 


different symmetiy adapted vectors out of ^ . It is to be noted that 

the choice of is completely arbitraiy and it is a matter of 
* 

convenience. By varying r,s and ^ , we arrive at vectors that have 
the symmetries of different irreducible representations. These are 

orthogonal, to one another and. with proper normalisation they constitute 

A 

the columns of the required V matrix ! which block diag-onolizes the 

matrices and g . This block diagonolization follows the properties 

of these column vectors and from the fact that the operators O' 

R 

commute with ^ and g . 

'••fe consider the specific case of the cubic group whose 

order N = 48 and whose character table is depicted in table 

o 

42 

III oelow in the notation of Hamermesh . Using equation (l III-3) 
for the simple cubic, body centred cubic face centred cubic and 

the polarization model, which are denoted by the subscripts S, B, F 
and P respectively, one gets with the aid of table III, 



s 


+ F. 
lu 

+ E 

g 

r 

B 

= 2 ilg 

+ F, 
lu 

+ F + A 

2g ^2u 

f 

TP 

= 2 A 

+ F 

+ E + F„ + Fo. 


r 

Ig 

lu 

g 2g 2u 

r 

p 

II 

® (2 

Ig lu 2g 2u ' 


where @ denotes the direct product of the irreducible 
repre sentati ons . 
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TABLE III 

CHjiB/CTER TABLE FOR 


Represen- 

tations. 

llota- „ 
tion®^^*^ 

E 

8 C 3 

3C| 

“2 

60, 

4 

I 



6 i 

5, 6S 
d 4 

A 

Ig 

s 

1 

1 

1 

1 

1 

1 

I 

1 

1 

1 

■^lu 

L 

1 

1 

1 

1 

1 . 

1 

-1 

-1 

-1 

-1 

^2g 

I 

1 

1 

1 

-1 

-1 

1 

1 

.1 

-1 

-1 

^2u 

^1 

1 

1 

1 

-1 

-1 

-1 

-1 

-1 

1 

1 

E 

g 


2 

-1 

2 

0 

0 

2 

-1 

2 

0 

0 


H 

2 

“1 

2 

0 

0 ‘ 

-2 

1 

-2 

0 

0 

F 

Ig 

G 

3 

0 

-1 

-1 

1 

3 

0 

-1 

-1 

1 

F 

^lu 

P 

3' 

0 

-1 

-1 

1 

-3 

0 

1 

1 

-1 

F 

2g 

^2 

t 

3 

0 

-1 

1 

-1 

3 

0 

-1 

1 

-1 

F 

?a 

1 

L__ 

3 

0 

-1 

1 

-1 

-3 

0 

1 

-1 

1 ^ I 



7 1 
9 3 
13 1 
27 0 


3 11 

111 
11 1 
1-3 0 


The matidcee ^ are worked out using the basis given in 

Hamerme suffer example. Now following the procedure described we arrive 

A 

, ^ at the following forms of V - matrices. in these, for each 

irreducible representation of dimensionality appearing once in the 
reduced form of P , one has columns.. If it appears a^ 

times then one has a^ x n^ columns. 
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and 


f \ 

The elements of 7-matrices given below are given in ternis oi 
numericel quantities given ty , 

b = (8)" ^ c = (24)" d = ^ , h = (12)'^, d' = (6)* 

h' = (2)" ^ 

(i) For the simple cubic lattice we get ^ 


A 

Yg = 


0 

0 

0 

0 

0 


Zk 


■Ig 


0 

0 

O' 

0 

0 

0 

d' 

h' 

0 

0 

h 

d 

d' 

0 

h' 

0 

h 

-d 

d' 

0 

0 

h» 

-2h 

0 

dl 

-h' 

0 

0 

h 

d 

d' 

0 

-h’ 

0 

h 

-d 

d' 

/ 

0 

0 

-h' 

-2h 

0 

J 





E 

g 


(4 III-4) 


Here the row and column indices, starting from upper left comer are 
(0,-0,0) , (1,0,0 ) , (0,1,0) , (0,0,1) , (-1,0,0) , (0,-l,0) and 
(0,0,-l) where the numbers refer to ( ^ >Di,n) values for the nearest 
neighbour sites according tc equation (3.1) of chapter IH. 

(ii) For the body centred cubic- lattice,- we have 



(AIII-5 ) 


with the row and column indices starting from theutper left comer 
as (0,0,0),(0,1,0)> (1,-1,0) , (0>-l,l) , (-1,1,-!), (0,-l,0) , (-1,1,0), 
(0,1, -l) and (l,-i,l) ; the numbers here referring to the ( L >ni,n) 
values for the nearest neighbour sites according to equations (4,1). 

(iii) B.nd for the face centred cubic lattice. 
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Vj. = 


/I 0 

j 0 h 
' 0 h 
0 h 
0 h 


0 h 


0 0 

b c 

b c 
0 -2c 

0 -2c 

-b c 
-b c 
b c 
b c 
0 -2c 

0 -2c 

-b c 
-b c 


6 0 0 

-d 0 0 

0 d 0 

0 -d 0 

0 0 d 

0 0 -d 


0 0 
0 0 


0 0 -d 


0 

0 

0 

b 

b 

0 

-b 

b 

0 

0 

-b 

b 

0 

-b 

-b 

b 

0 

-b 

-b 

0 

-b 

-b 

-b 

0 

b 

-b 

0 

0 

b 

-b 

0 

b 

b 

-b 

0 

b 

b 

0 

b 


iUin-6) 


with the row and column indices starting from the upper left comer 

as (0,0,0) , (1,0,0) , (C,l,-1) , (0,1,0) , (l,0,-l) , (-1,1,0) ,(0,0,l) 

(-1,0,0) , (0,-l,l) , (0,-1, 0), (-1,0,1) , (1,-1,0) and (0,0,-l) . 

A 

Finally the matrix Vp for the model with polarization is a 
27 X 27 matrix which can be partitioned into 9x9 submatiices as shown 
in equation (5.53) of chapter V. 

n /V A 


Vp = 


"'ll ^12 ^13 

A /V A 

''^21 ^22 ^23 

'V A /\ 

V V V 

31 32 33 


(AIII-7) 


Each of these submatrices have row and column indices the same as 
A 

described for Vp above ( equation (aIII- 5) ). The older of columns 
A 

of Vp will coriespond to the following order of the irreducible 
representations of Oj^ appearing in the reduced fom of . 

lu ' Ig ' g ' 2g ' 2u ’ 2u ' 10 u 

This leads to the order of blocks described in chapter V. 

A 

Explicitly, the matrices V^^ ,(i,j = 1,2,3) of equation 
( AIII-7) can be lArritten ty using the follovring nine dimensional 


vectors . 




Xi = (fi)A 


Z. =l><i 

1 


'f-± 


for i = 3,4, 


B^r writing these vectors side by side with appropriate factors and 

■ ^ 

signs , the subcatricos’ V- . .are. 


V = 

(Xi 

■^0 

x^ 

X 

X 

X 

X 

X 

X 

) 

11 


2 

0 

0 

0 

8 

0 

0 

7 


Y = 
12 

(h 

^4 

24 

^3 

^0 



^0 

^9 

) 

^3 "" 

(^6 

^7 


"8 

"0 

"5 

^3 

^6 

^6 

) 

'^21 = 





^2 


^0 

^0 

^6 

) 

^•"22 ^ 

(b 

-^3 

23 

^4 

^0 


^9 

^5 

^0 

) 

V = 


-X 

-X 

X 

-X 

X 

-X 

z 

-X 

) 

23 

7 

6 

8 

0 

5 

0 

4 

7 

7 


^'31 = 


^0 


"0 

^0 

^6 




) 

^32 = 

(Y 

'■5 

"0 

-^25 

^0 

^9 


^0 

Xj 

^0 

) 

Y = 


X 


-X 

X 

-X 

X 

-22 

X 

) 

33 

8 

0 

7 

6 

3 

4 

0 

8 

0 



^Aflien all these are put together in equation (aIII- 7) one gets the 

A 

full 27 X 27 matrix Yp which is used in chapter V. 
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